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The pure-gravity sector of the minimal Standard-Model Extension is studied in the limit of Rie- 
mann spacetime. A method is developed to extract the modified Einstein field equations in the limit 
of small metric fluctuations about the Minkowski vacuum, while allowing for the dynamics of the 
20 independent coefficients for Lorentz violation. The linearized efi^ective equations are solved to 
obtain the post-newtonian metric. The corresponding post-newtonian behavior of a perfect fluid is 
studied and applied to the gravitating many-body system. Illustrative examples of the methodology 
are provided using bumblebee models. The implications of the general theoretical results are stud- 
ied for a variety of existing and proposed gravitational experiments, including lunar and satellite 
laser ranging, laboratory experiments with gravimeters and torsion pendula, measurements of the 
spin precession of orbiting gyroscopes, timing studies of signals from binary pulsars, and the classic 
tests involving the perihelion precession and the time delay of light. For each type of experiment 
considered, estimates of the attainable sensitivities are provided. Numerous effects of local Lorentz 
violation can be studied in existing or near-future experiments at sensitivities ranging from parts in 
10* down to parts in lO^'^. 
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I. INTRODUCTION 

At the classical level, gravitational phenomena are well 
described by general relativity, which has now survived 
nine decades of experimental and theoretical scrutiny. 
In the quantum domain, the Standard Model of parti- 
cle physics offers an accurate description of matter and 
nongravitational forces. These two field theories pro- 
vide a comprehensive and successful description of na- 
ture. However, it remains an elusive challenge to find a 
consistent quantum theory of gravity that would merge 
them into a single underlying unified theory at the Planck 
scale. 

Since direct measurements at the Planck scale are in- 
feasible, experimental clues about this underlying theory 
are scant. One practical approach is to search for prop- 
erties of the underlying theory that could be manifest 
as suppressed new physics effects, detectable in sensi- 
tive experiments at attainable energy scales. Promising 
candidate signals of this type include ones arising from 
minuscule violations of Lorentz symmetry [1, 2]. 

Effective field theory is a useful tool for describing ob- 
servable signals of Lorentz violation [3]. Any realistic 
effective field theory must contain the Lagrange densi- 
ties for both general relativity and the Standard Model, 
possibly along with suppressed operators of higher mass 
dimension. Adding also all terms that involve operators 
for Lorentz violation and that are scalars under coordi- 
nate transformations results in an effective field theory 
called the Standard-Model Extension (SME). The lead- 
ing terms in this theory include those of general relativ- 
ity and of the minimally coupled Standard Model, along 
with possible Lorentz-violating terms constructed from 
gravitational and Standard-Model fields. 

Since the SME is founded on well-established physics 
and constructed from general operators, it offers an ap- 
proach to describing Lorentz violation that is largely in- 
dependent of the underlying theory. Experimental pre- 



dictions of realistic theories involving relativity modifi- 
cations are therefore expressible in terms of the SME by 
specifying the SME coefficient values, fn fact, the explicit 
form of all dominant Lorentz-violating terms in the SME 
is known [4]. These terms consist of Lorentz-violating 
operators of mass dimension three or four, coupled to 
coefficients with Lorentz indices controlling the degree 
of Lorentz violation. The subset of the theory contain- 
ing these dominant Lorentz-violating terms is called the 
minimal SME. 

Since Lorentz symmetry underlies both general rela- 
tivity and the Standard Model, experimental searches for 
violations can take advantage either of gravitational or of 
nongravitational forces, or of both. In the present work, 
we initiate an SME-based study of gravitational experi- 
ments searching for violations of local Lorentz invariance. 
To restrict the scope of the work to a reasonable size 
while maintaining a good degree of generality, we limit 
attention here to the pure-gravity sector of the minimal 
SME in Riemann spacetime. This neglects possible com- 
plexities associated with matter-sector effects and with 
Riemann-Cartan [5] or other generalized spacetimes, but 
it includes all dominant Lorentz-violating signals in ef- 
fective action-based metric theories of gravity. 

The Minkowski-spacetime limit of the minimal SME 
[6] has been the focus of various experimental studies, 
including ones with photons [7-9], electrons [10-12], pro- 
tons and neutrons [13-15], mesons [16], muons [17], neu- 
trinos [18], and the Higgs [19]. To date, no compelling 
evidence for nonzero coefficients for Lorentz violation has 
been found, but only about a third of the possible sig- 
nals involving light and ordinary matter (electrons, pro- 
tons, and neutrons) have been experimentally investi- 
gated, while some of the other sectors remain virtually 
unexplored. Our goal here is to provide the theoretical 
basis required to extend the experimental studies into the 
post-newtonian regime, where asymptotically Minkowski 
spacetime replaces the special-relativistic limit. 
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Nonzero SME coefficients for Lorentz violation could 
arise via several mechanisms. It is convenient to dis- 
tinguish two possibilities, spontaneous Lorentz violation 
and explicit Lorentz violation. If the Lorentz violation 
is spontaneous [20], the SME coefficients arise from un- 
derlying dynamics, and so they must be regarded as 
fields contributing to the dynamics through the varia- 
tion of the action. In contrast, if the Lorentz viola- 
tion is explicit, the SME coefficients for Lorentz vio- 
lation originate as prescribed spacetime functions that 
play no dynamical role. Remarkably, the geometry of 
Riemann-Cartan spacctimes, including the usual Rie- 
mann limit, is inconsistent with explicit Lorentz viola- 
tion [4]. In principle, a more general non-Riemann ge- 
ometry such as a Finsler geometry might allow for ex- 
plicit violation [4, 21], but this possibility remains an 
open issue at present. We therefore limit attention to 
spontaneous Lorentz violation in this work. Various sce- 
narios for the underlying theory are compatible with a 
spontaneous origin for Lorentz violation, including ones 
based on string theory [20], noncommutative field theo- 
ries [22], spacetime- varying fields [23], quantum gravity 
[24], random-dynamics models [25], multiverses [26], and 
brane- world scenarios [27]. 

Within the assumption of spontaneous Lorentz break- 
ing, the primary challenge to extracting the post- 
newtonian physics of the pure-gravity minimal SME lies 
in accounting correctly for the fluctuations around the 
vacuum values of the coefficients for Lorentz violation, in- 
cluding the massless Nambu-Goldstone modes [28]. Ad- 
dressing this challenge is the subject of Sec. II of this 
work. The theoretical essentials for the analysis are pre- 
sented in Sec. II A, while Sec. II B describes our method- 
ology for obtaining the effective Unearized field equations 
for the metric fluctuations in a general scenario. The 
post-newtonian metric of the pure-gravity minimal SME 
is obtained in Sec. Ill A, and it is used to discuss modifi- 
cations to perfect-fiuid and many-body dynamics in Sec. 
IIIB. In recent decades, a substantial effort has been in- 
vested in analysing weak-field tests of general relativity 
in the context of post-newtonian expansions of an arbi- 
trary metric, following the pioneering theoretical works of 
Nordtvedt and Will [29, 30]. Some standard and widely 
used forms of these expansions are compared and con- 
trasted to our results in Sec. IIIC. The theoretical part 
of this work concludes in Sec. IV, where the key ideas for 
our general methodology are illustrated in the context of 
a class of bumblebee models. 

The bulk of the present paper concerns the implica- 
tions for gravitational experiments of the post-newtonian 
metric for the pure-gravity sector of the minimal SME. 
This issue is addressed in Sec. V. To keep the scope 
reasonable, we limit attention to leading and subleading 
effects. We identify experiments that have the potential 
to measure SME coefficients for Lorentz violation, and 
we provide estimates of the attainable sensitivities. The 
analysis begins in Sec. VA with a description of some 
general considerations that apply across a variety of ex- 



periments, along with a discussion of existing bounds. In 
Sec. VB, wc focus on experiments involving laser rang- 
ing, including ranging to the Moon and to artificial satel- 
lites. Section V C studies some promising laboratory ex- 
periments on the Earth, including gravimeter measure- 
ments of vertical acceleration and torsion-pendulum tests 
of horizontal accelerations. The subject of Sec. V D is the 
precession of an orbiting gyroscope, while signals from 
binary pulsars are investigated in Sec. VE. The role of 
the classic tests of general relativity is discussed in Sec. 
VP. To conclude the paper, a summary of the main 
results is provided in Sec. VI, along with a tabulation 
of the estimated attainable experimental sensitivities for 
the SME coefficients for Lorentz violation. Some details 
of the orbital analysis required for our considerations of 
laser ranging are relegated to Appendix A. Throughout 
this work, we adopt the notation and conventions of Ref. 
[4]. 



II. THEORY 
A. Basics 

The SME action with gravitational couplings is pre- 
sented in Ref. [4]. In the general case, the geomet- 
ric framework assumed is a Riemann-Cartan spacetime, 
which allows for nonzero torsion. The pure gravitational 
part of the SME Lagrange density in Riemann-Cartan 
spacetime can be viewed as the sum of two pieces, one 
Lorentz invariant and the other Lorentz violating: 

■^gravity = ^LI + -^LV- (1) 

All terms in this Lagrange density are invariant un- 
der observer transformations, in which all fields and 
backgrounds transform. These include observer local 
Lorentz transformations and observer diffeomorphisms 
or general coordinate transformations. The piece £li 
also remains invariant under particle transformations, in 
which the localized fields and particles transform but the 
backgrounds remain fixed. These include particle local 
Lorentz transformations and particle diffeomorphisms. 
However, for vanishing fiuctuations of the coefficients for 
Lorentz violation, the piece £lv changes under particle 
transformations and thereby breaks Lorentz invariance. 

The Lorentz-invariant piece Chi is a series in pow- 
ers of the curvature R^Xfiu, the torsion T\^i,, the co- 
variant derivatives D^,, and possibly other dynamical 
fields determining the pure-gravity properties of the the- 
ory. The leading terms in Cia are usually taken as 
the Einstein-Hilbert and cosmological terms in Riemann- 
Cartan spacetime. The Lorentz-violating piece £lv is 
constructed by combining coefficients for Lorentz viola- 
tion with gravitational field operators to produce indi- 
vidual terms that are observer invariant under both local 
Lorentz and general coordinate transformations. The ex- 
plicit form of this second piece can also be written as a 
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series in the curvature, torsion, covariant derivative, and 
possibly other fields: 

+e{kDTT^^''D^n^, + . . . , (2) 

where e is the determinant of the vierbein e,,^. The 
coefficients for Lorentz violation (kT)^^" , {knY^^^ , 
{kTq^)"^^''^'^''\ {kD-rY^^'^ can vary with spacctime posi- 
tion. Since particle local Lorentz violation is always ac- 
companied by particle diffeomorphism violation [28] , the 
coefficients for Lorentz violation also control diffeomor- 
phism violation in the theory. 

In the present work, we focus on the Riemann- 
spacetime limit of the SME, so the torsion is taken to 
vanish. We suppose that the Lorentz-invariant piece of 
the theory is the Einstcin-Hilbcrt action, and we also 
restrict attention to the leading-order Lorentz-violating 
terms. The gravitational terms that remain in this limit 
form part of the minimal SME. The basic features of the 
resulting theory are discussed in Ref. [4], and those rele- 
vant for our purposes are summarized in this subsection. 

The effective action of the minimal SME in this limit 
can be written as 

5 = 5eh + ^LV + S'. (3) 

The first term in (3) is the Einstein-Hilbert action of 
general relativity. It is given by 

S^n = Y^jd^xe{R-2K), (4) 

where R is the Ricci scalar, A is the cosmological con- 
stant, and K = SttG. As usual, in the present context of a 
Riemann spacetime, the independent degrees of freedom 
of the gravitational field are contained in the metric g^^. 
Since we are ultimately focusing on the post-newtonian 
limit of (3), in which the effects of A are known to be 
negligible, we set A = for the remainder of this work. 

The second term in Eq. (3) contains the leading 
Lorentz-violating gravitational couplings. They can be 
written as 

5lv = ^ y d^xe{-uR + s^'^i?^, + 1''^^^ C ^Xf..) ■ (5) 

In this equation, i?^^^ is the trace-free Ricci tensor and 
CkXhv is the Weyl conformal tensor. The coefficients for 
Lorentz violation s^" and t'^'^^^ inherit the symmetries 
of the Ricci tensor and the Riemann curvature tensor, 
respectively. The structure of Eq. (5) implies that s^" 
can be taken traceless and that the various traces of t"-^i^'^ 
can all be taken to vanish. It follows that Eq. (5) contains 
20 independent coefficients, of which one is in u, 9 are in 
si^", and 10 are in f^^i"" . 

The coefficients w, s'"^, and t'^^^^'^ typically depend on 
spacetime position. Their nature depends in part on the 



origin of the Lorentz violation. As mentioned in the 
introduction, explicit Lorentz violation is incompatible 
with Riemann spacetime [4] . We therefore limit attention 
in this work to spontaneous Lorentz violation in Riemann 
spacetime, for which the coefficients u, s'"', t'^^i^" are dy- 
namical fields. Note that spontaneous local Lorentz vi- 
olation is accompanied by spontaneous diffeomorphism 
violation, so as many as 10 symmetry generators can be 
broken through the dynamics, with a variety of interest- 
ing attendant phenomena [28]. Note also that u, s'"^, 
^kXixv jj^g^y composites of fields in the underlying the- 
ory. Examples for this situation are discussed in Sec. IV. 

The third term in Eq. (3) is the general matter action 
S' . In addition to determining the dynamics of ordinary 
matter, it includes contributions from the coefficients u, 
s'"', t'^^^'^ , which for our purposes must be considered in 
some detail. The action S' could also be taken to include 
the SME terms describing Lorentz violation in the matter 
sector. These terms, given in Ref. [4], include Lorentz- 
violating matter-gravity couplings with potentially ob- 
servable consequences, but addressing these effects lies 
beyond the scope of the present work. Here, we focus 
instead on effects from the gravitational and matter cou- 
plings of the coefficients u, s^", t'^^t^'^ in Eq. (5). 

Variation with respect to the metric g^^ while holding 
u, s'^^ , and t'^^^^'^ fixed yields the field equations 

Qixu _ ^r^Rstuyu ^ KiTgY". (6) 
In this expression, 

-D^Df:)t''''''^ - DaDpt'"^'"^, (7) 

while the general matter energy-momentum tensor is de- 
fined as usual by 

^e{T,r^6C'/5g^., (8) 

where C is the Lagrange density of the general matter 
action S'. 



B. Linearization 

One of the central goals of this work is to use the SME 
to obtain the newtonian and leading post-newtonian cor- 
rections to general relativity induced by Lorentz viola- 
tions. For this purpose, it suffices to work at linear order 
in metric fluctuations about a Minkowski background. 
We can therefore adopt the usual asymptotically inertial 
coordinates and write 

9^,u = riiiu + hij,^. (9) 

In this subsection, we derive the effective linearized field 
equations for h^i, in the presence of Lorentz violation. 
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1. Primary linearization 

A key issue is the treatment of the dynamics of the 

coefficient fields u, s'"', and t'^^t^'^ . As described above, 
these are assumed to induce spontaneous violation of lo- 
cal Lorentz invariance and thereby acquire vacuum ex- 



pectation values u, s^^ , and t 



respectively. Denot- 



ing the field fluctuations about these vacuum solutions 
as M, s'^" , and t'^^i^" . we can write 



u 



u + u, 



+ s 



t 



(10) 



Note that the fluctuations include as a subset the Nambu- 
Goldstone modes for local Lorentz and diffeomorphism 
violation, which are described in Ref. [28]. For present 
purposes it suffices to work at linear order in the fluc- 
tuations, so in what follows nonlinear terms at 
0{hu), O(s^), etc., are disregarded, and we adopt the 
standard practice of raising and lowering indices on lin- 
ear quantities with the Minkowski metric t^^j/. 

Deriving the effective linearized field equations for h^i, 
involves applying the expressions (9)-(10) in the asymp- 
totically inertial frame. It also requires developing meth- 
ods to account for effects on h^^i, due to the fluctuations 
u, 8^^" , i'^^t^'^. In fact, as we show below, five key assump- 
tions about the properties of these fluctuations suffice for 
this purpose. 

The first assumption concerns the vacuum expectation 
values. We assume (i): the vacuum values u, s'^" , t'^'^^" 
are constant in asymptotically inertial cartesian coordi- 
nates. Explicitly, we take 



daU = 



dryt 



■KXfll/ 



0, 
0, 
0. 



(11) 



More general conditions could be adopted, but assump- 
tion (i) ensures that translation invariance and hence 
energy-momentum conservation are preserved in the 
asymptotically Minkowski regime. The reader is cau- 
tioned that this assumption is typically different from 
the requirement of covariant constancy. 

The second assumption ensures small Lorentz- 
violating effects. We suppose (ii): the dominant effects 
are linear in the vacuum values u, s^", J'^^^'^ ^ This as- 
sumption has been widely adopted in Lorentz-violation 
phenomenology. The basic reasoning is that any Lorentz 
violation in nature must be small, and hence lineariza- 
tion typically suffices. Note, however, that in the present 
context this assumption applies only to the vacuum val- 
ues of the coefficients for Lorentz violation in the SME 
action (3). Since these coefficients arc related via un- 
determined coupling constants to the vacuum values of 
the dynamical fields in the underlying theory, assump- 
tion (ii) provides no direct information about the sizes of 
the latter. 



With these first two assumptions, we can extract the 
linearized version of the field equations (7) in terms of h^^ 
and the fluctuations m, s^^^ , and t'^^i^" . Some calculation 
shows that the linearized trace-reversed equations can be 
written in the form 



(12) 



In this expression, {Sg)^^ is the trace-reversed energy- 
momentum tensor deflned by 



a 1 



(13) 



where (Tg)^^ is the linearized energy-momentum tensor 
containing the energy-momentum density of both con- 
ventional matter and the fluctuations u, s^^ , and t'^^^^^ . 
Also, the terms A^iu and Bfj,^ in Eq. (12) are given by 



2 D-S^iy 



(14) 



and 



(15) 



The connection coefficients appearing in Eq. (14) are lin- 
earized Christoffel symbols, where indices are lowered 
with the Minkowski metric as needed. The terms R, Ra/s, 
and Ra^-yS appearing in Eqs. (12)-(15) and elsewhere be- 
low are understood to be the linearized Ricci scalar. Ricci 
tensor, and Riemann curvature tensor, respectively. 



2. Treatment of the energy-momentum tensor 



To generate the effective equation of motion for /i^j^ 
alone, the contributions from the fluctuations u, s'^", 
^kXixv jjjyg^ expressed in terms of hnv, its derivatives, 
and the vacimm values u, s^'^ , ^'^^a'^ general, this is a 
challenging task. We adopt here a third assumption that 
simplifles the treatment of the dynamics sufflciently to 
permit a solution while keeping most interesting cases. 
We assume (iii): the fluctuations u. s^"\ f^^i^'^ have no 
relevant couplings to conventional m,atter. This assump- 
tion is standard in alternate theories of gravity, where it 
is desired to introduce new fields to modify gravity while 
maintaining suitable matter properties. It follows that 
{Sg)^i, can be split as 



{Sg)fj,v — {,Sm)ixv + {,Sstu)nvi 



(16) 



where (S'm)/^i> is the trace-reversed energy-momentum 
tensor for conventional matter. 
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To gain insight about assumption (iii), consider the 
contributions of the fluctuations w,, S^"^, f'^V'' ^j^g 
trace-reversed energy-momentum tensor {Sg)i^i,. A pri- 
ori, the fields u, s^" , t'^^t^" can have couphngs to mat- 
ter currents in the action S", which would afi'ect the 
energy-momentum contribution of u, s^^ , fK\fj.i^ _ How- 
ever, through the vacuum values u, s^'^ , J'^^^'^ ^ these 
couplings would also induce coefficients for Lorentz vi- 
olation in the SME matter sector, which are generically 
known to be small from nongravitational experiments. 
The couplings of the fluctuations u, s^^ , ^^A/ii/ ^j^g 
matter currents are therefore also generically expected to 
be small, and so it is reasonable to take the the matter 
as decoupled from the fluctuations for present purposes. 

Possible exceptions to the decoupling of the fluctua- 
tions and matter can arise. For example, in the case of 
certain bumblebee models, including those described in 
Sec. IV, the fluctuations s^'^ can be written in terms of a 
vacuum value and a propagating vector field. This vec- 
tor field can be identified with the photon in an axial 
gauge [28]. A sufficiently large charge current could 
then generate significant fluctuations s'*'', thereby com- 
peting with gravitational effects. However, for the var- 
ious systems considered in this work, the gravitational 
fields dominate over the electrodynamic fields by a con- 
siderable amount, so it is again reasonable to adopt as- 
sumption (iii). 

Given the decomposition (16), the task of decoupling 
the fluctuations u, s^^^ , ^kXhu (S'g)^,^ reduces to ex- 
pressing the partial energy-momentum tensor {Sstu)^iv 
in terms of h^v, its derivatives, and the vacuum values 
u, s'"^, t'^^^'^ . For this purpose, we can apply a set of 
four identities, derived from the traced Bianchi identi- 
ties, that are always satisfied by the gravitational energy- 
momentum tensor [4]. The linearized versions of these 
conditions suffice here. They read 

^^^^{Tg)^.. = -s^^^^i?.. - 2r^^^ dsRafiju- (17) 

Using the fact that {TM)fj.u is separately conserved, one 
can show that the four conditions (17) are satisfied by 

^^{S stu) 1X1/ = ~2s"^^i?j^')„ + ^S^i,R + l^i^iyS"^ Raij 

+4r/,i?„/3 + S^,. (18) 

In this equation, the term E^,^ obeys 

d^{Y.^, - i77^.E«J = 0. (19) 

It represents an independently conserved piece of the lin- 
earized energy-momentum tensor that is undetermined 
by the conditions (17). 

For calculational purposes, it is convenient in what fol- 
lows to adopt assumption (iv): the independently con- 
served piece of the trace-reversed energy momentum ten- 
sor vanishes, S^,^ = 0. Since E^^ is independent of other 
contributions to {Sstu)nv, one might perhaps suspect that 



it vanishes in most theories, at least to linear order. How- 
ever, theories with nonzero E^j^ do exist and may even be 
generic. Some simple examples are discussed in Sec. IV. 
Nonetheless, it turns out that assumption (iv) suffices for 
meaningful progress because in many models the nonzero 
contributions from E^^^ merely act to scale the effective 
linearized Einstein equations relative to those obtained in 
the zero-E^jy limit. Models of this type are said to vio- 
late assumption (iv) weakly, and their linearized behavior 
is closely related to that of the zero-E^i, limit. In con- 
trast, a different behavior is exhibited by certain models 
with ghost kinetic terms for the basic fields. These have 
nonzero contributions to E^^ that qualitatively change 
the behavior of the effective linearized Einstein equations 
relative to the zero-E^i^j, limit. Models with this feature 
are said to violate assumption (iv) strongly. It is reason- 
able to conjecture that for propagating modes a strongly 
nonzero E^,^ is associated with ghost fields, but estab- 
lishing this remains an open issue lying outside the scope 
of the present work. 



3. Decoupling of fluctuations 

At this stage, the trace-reversed energy-momentum 

tensor {Sg)fj_i, in Eq. (12) has been linearized and ex- 
pressed in terms of h^j^, its derivatives, and the vacuum 
values u, 's'^" , j^x^v ^^j^^ term B^j/ in Eq. (12) already 
has the desired form, so it remains to determine Afj,v 
The latter explicitly contains the fluctuations s^", t'^^t^" ^ 
and u. By assumption (iii), the fluctuations couple only 
to gravity, so it is possible in principle to obtain them 
as functions of h^^, alone by solving their equations of 
motion. It follows that A^u can be expressed in terms 
of derivatives of these functions. Since the leading-order 
dynamics is controlled by second-order derivatives, the 
leading-order result for A^,,^ is also expected to be scicond 
order in derivatives. We therefore adopt assumption (v): 
the undetermined terms in Afiv are constructed from lin- 
ear combinations of two partial spacetime derivatives of 
h^i, and the vacuum values ri^i,, u, 's^'^ , andt!^^^^. More 
explicitly, the undetermined terms take the generic form 
^^iij"^'^^ dadph^s- This assumption ensures a smooth 
match to conventional general relativity in the limit of 
vanishing coefficients M^^"^'>"'. 

To constrain the form of Af_tu, we combine assiimp- 
tion (v) with invariance properties of the action, notably 
those of diffeomorphism symmetry. Since we are consid- 
ering spontaneous symmetry breaking, which maintains 
the symmetry of the full equations of motion, the origi- 
nal particle diffeomorphism invariance can be applied. In 
particular, it turns out that the linearized particle diffeo- 
morphism transformations leaves the linearized equations 
of motion (12) invariant. To proceed, note first that the 
vacuum values u, s^"^ , J'^^^^'^ and rj^,^ are invariant under 
a particle diffeomorphism parametrized by [28] , while 
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the metric fluctuation transforms as 

h^u hf,^ - dfj,^^ - d^^fj,. (20) 

The form of Eq. (18) then implies that {Sg)^^ is invariant 
at linear order, and from Eq. (15) we find B^v is also 
invariant. It then follows from Eq. (12) that Afj,v must 
be particle diffeomorphism invariant also. 

The invariance of Afiu can also be checked directly. 
Under a particle diffeomorphism, the induced transfor- 
mations on the fluctuations s'"', t'^^t^'^^ u can be derived 
from the original transformation of the fields u, s^^ , and 
^kX^iv ^Yic same manner that the induced transforma- 
tion of is obtained from the transformation for g^n. 
We find 

U ^ U. (21) 

Using these transformations, one can verify explicitly 
that A^v is invariant at leading order under particle dif- 
feomorphisms. 

The combination of assumption (v) and the imposition 
of particle diffeomorphism invariance suffices to extract 
a covariant form for A^j,ij involving only the metric fluc- 
tuation h^j^. After some calculation, we find A^u takes 
the form 

A^v = —2auRfj,^ + S°'^Ra(^^)i3 — s^^^^R^^a 

-hT'^^ (^^R^^^ap- (22) 

Some freedom still remains in the structure of A^v, as 
evidenced in Eq. (22) through the presence of arbitrary 
scaling factors a and b. These can take different values 
in distinct explicit theories. 

We note in passing that, although the above derivation 
makes use of particle diffeomorphism invariance, an alter- 
native possibility exists. One can instead apply observer 
diffeomorphism invariance, which is equivalent to invari- 
ance under general coordinate transformations. This in- 
variance is unaffected by any stage of the spontaneous 
symmetry breaking, so it can be adapted to a version of 
the above reasoning. Some care is required in this proce- 
dure. For example, under observer diffeomorphisms the 
vacuum values u, s^^ , ji^^t^^ transform nontrivially, and 
the effects of this transformation on Eqs. (11) must be 
taken into accoimt. In contrast, the fluctuations .1^'', 
^KXfii^^ ~ unchanged at leading order under observer 
diffeomorphisms. In any case, the result (22) provides 
the necessary structure of A^v when the assumptions (i) 
to (v) are adopted. 



4- Effective linearized Einstein equations 

The flnal effective Einstein equations for the metric 

fluctuation h^^ are obtained upon inserting Eqs. (15), 
(18), and (22) into Eq. (12). We can arrange the equa- 
tions in the form 

Ri.. = k{Sm)>.. + -^% + ^1. + <pI,.. (23) 

The quantities on the right-hand side are given by 

= uRf,^, 

^Itu = 2i'*'^'''(^i?i.)^a/3 + 2t^ fRaP 

= 0. (24) 

Each of these quantities is independently conserved, as 

required by the linearized Bianchi identities. Each is also 
invariant under particle diffeomorphisms, as can be ver- 
ified by direct calculation. The vanishing of at the 
linearized level is a consequence of these conditions and of 
the index structure of the coefficient t"^^^ ^ which implies 
the identity [31] t^^^'^^ 6''^^^ = 0. 

In Eq. (23), the coefficients of <i>*,^, are taken 
to be unity by convention. In fact, scalings can arise from 
the terms in Eq. (22) or, in the case of models weakly 
violating assumption (iv), from a nonzero S^^/ term in 
Eq. (18). However, these scalings can always be absorbed 
into the definitions of the vacuum values u, s^" , f'^^'^'^ _ 
In writing Eq. (23) we have implemented this rescaling of 
vaciium values, since it is convenient for the calculations 
to follow. However, the reader is warned that as a result 
the vacuum values of the fields u, s^^ , t'^^^^ appearing in 
Eq. (5) differ from the vacuum values in Eqs. (24) by a 
possible scaling that varies with the specific theory being 
considered. 

Since Eq. (23) is linearized both in and in the 
vacuum values, the solution for h^v can be split into the 
sum of two pieces, 

V = C + V' (25) 

one conventional and one depending on the vacuum val- 
ues. The first term h^^ is defined by the requirement 
that it satisfy 

= h{Sm),iu; (26) 

which are the standard linearized Einstein equations of 
general relativity. The second piece controls the devia- 
tions due to Lorentz violation. Denoting by R^i, the Ricci 
tensor constructed with hnv, it follows that this second 
piece is determined by the expression 

^/^^ = + (27) 
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where it is understood that the terms on the right-hand 
side are those of Eqs. (24) in the hmit /ij^^ 



so 



that only terms at linear order in Lorentz violation are 
kept. 

Equation (27) is the desired end product of the lin- 
earization process. It determines the leading corrections 
to general relativity arising from Lorentz violation in a 
broad class of theories. This includes any modified the- 
ory of gravity that has an action with leading-order con- 
tributions expressible in the form (5) and satisfying as- 
sumptions (i)-(v). Note that the fields u, s^^" , i^Xnv 
be composite, as occurs in the bumblebee examples dis- 
cussed in Sec. IV. Understanding the implications of Eq. 
(27) for the post-newtonian metric and for gravitational 
experiments is the focus of the remainder of this work. 



III. POST-NEWTONIAN EXPANSION 

This section performs a post-newtonian analysis of the 
linearized effective Einstein equations (27) for the pure- 
gravity sector of the minimal SME. We first present the 
post-newtonian metric that solves the equations. Next, 
the cqiiations of motion for a perfect fluid in this met- 
ric are obtained. Applying them to a system of massive 
self-gravitating bodies yields the leading-order accelera- 
tion and the lagrangian in the point-mass limit. Finally, 
a comparison of the post-newtonian metric for the SME 
with some other known post-newtonian metrics is pro- 
vided. 



A. Metric 

Following standard techniques [32] , we expand the lin- 
earized effective Einstein equations (26) and leading- 
order corrections (27) in a post-newtonian series. The 
relevant expansion parameter is the typical small velocity 
V oi a body within the dynamical system, which is taken 
to be 0(1). The dominant contribution to the metric 
fluctuation /i^y is the newtonian gravitational potential 
U. It is second order, 0{2) = ?a GM /r, where M is 
the typical body mass and f is the typical system dis- 
tance. The source of the gravitational field is taken to 
be a perfect fiuid, and its energy-momentum tensor is 
also expanded in a post-newtonian series. The dominant 
term is the mass density p. The expansion for h^^ begins 
at 0(2) because the leading-order gravitational equation 
is the Poisson equation V^?7 = —AnGp. 

The focus of the present work is the dominant Lorentz- 
violating effects. We therefore restrict attention to the 
newtonian 0(2) and sub-newtonian 0(3) corrections in- 
duced by Lorentz violation. For certain experimental 
applications, the 0(4) metric fiuctuation hoo would in 
principle be of interest, but deriving it requires solving 
the sub- linearized theory of Sec. II A and lies beyond the 
scope of the present work. 



As might be expected from the form of in Eq. (24), 
which involves a factor u multiplying the Ricci tensor, a 
nonzero u acts merely to scale the post-newtonian metric 
derived below. Moreover, since the vacuum value u is a 
scalar under particle Lorentz transformations and is also 
constant in asymptotically inertial coordinates, it plays 
no direct role in considerations of Lorentz violation. For 
convenience and simplicity, we therefore set u = in 
what follows. However, no assumptions are made about 
the sizes of the coefficients for Lorentz violation, other 
than assuming they are sufficiently small to validate the 
perturbation techniques adopted in Sec. II B. In terms 
of the post-newtonian bookkeeping, we treat the coeffi- 
cients for Lorentz violation as 0(0). This ensures that we 
keep all possible Lorentz- violating corrections implied by 
the linearized field equations (23) at each post-newtonian 
order considered. 

The choice of observer frame of reference affects the 
coefficients for Lorentz violation and must therefore be 
specified in discussions of physical effects. For immediate 
purposes, it suffices to assume that the reference frame 
chosen for the analysis is approximately asymptotically 
inertial on the time scales relevant for any experiments. 
In practice, this implies adopting a reference frame that 
is comoving with respect to the dynamical system un- 
der consideration. The issue of specifying the observer 
frame of reference is revisited in more detail as needed in 
subsequent sections. 

As usual, the development of the post-newtonian se- 
ries for the metric involves the introduction of certain 
potentials for the perfect fluid [30]. For the pure-gravity 
sector of the minimal SME taken at 0(3), we find that 
the following potentials are required: 



U 
Ijjk 



,p {x',t) 
R ' 



G j d^x- 
G j d^x'^-^^^p^, 



= G j d^x' 



G 



G j d^x' 
? j d^x- 



,p {x',t)v^{x',t) 
R 

,3 ,p{x' ,t)vHS',t)R''R^ 
R^ 

,p{x \t)v^{x\t)R''R^ 



G I ,s,>Pi^,t)v-ix',m-R^R^R^^ ^2^^ 



i?5 



x'^ and R ■ 



with the euclidean 



where R^ = 
norm. 

The potential U is the usual newtonian gravitational 
potential. In typical gauges, the potentials , oc- 
cur in the post-newtonian expansion of general relativ- 
ity, where they control various gravitomagnetic effects. 
In these gauges, the potentials U^'', X^^\ and F-''"' lie 
beyond general relativity. To our knowledge the latter 
two, X^^^ and Y^^\ have not previously been considered 
in the literature. However, they are needed to construct 
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the contributions to the 0(3) metric arising from general 
leading-order Lorentz violation. 
A 'superpotential' x defined by 

X = -G j (fx'p{x',t)R (29) 

is also used in the literature [30]. ft obeys the identities 

dodjx = -W\ (30) 

In the present context, it is convenient to introduce two 
additional superpotentials. We define 

= G J d^x'p{x',ty{x',t)R (31) 

and 

= G j (fx'p{o^,ty{o^,t)R>R. (32) 

These obey several useful identities including, for exam- 
ple, 

djdkX^ = P'^V^ - 

(33) 

In the latter equation, the parentheses denote total sym- 
metrization with a factor of 1/3. 

In presenting the post-newtonian metric, it is neces- 
sary to fix the gauge. In our context, it turns out that 
calculations can be substantially simplified by imposing 
the following gauge conditions: 

9j9jk = \dk{9jj - 9oo)- (34) 

It is understood that these conditions apply to 0(3). Al- 
though the conditions (34) appear superficially similar to 
those of the standard harmonic gauge [32], the reader is 
warned that in fact they differ at 0(3). 

With these considerations in place, direct calculation 
now yields the post-newtonian metric at 0(3) in the cho- 
sen gauge. The procedure is to break the effective lin- 
earized equations (27) into temporal and spatial compo- 
nents, and then to use the usual Einstein equations to 
eliminate the pieces h^^ of the metric on the right-hand 
side in favor of the potentials (28) in the chosen gauge, 
keeping appropriate track of the post-newtonian orders. 
The resulting second-order differential equations for /i^,^ 
can be solved in terms of the potentials (27) . 

After some work, we find that the metric satisfying 



Eqs. (26) and (27) can be written at this order as 

£,00 = -1 + 2U + 2,-^°U + a>''W^ - ^-^^V^ +0{A), 

(35) 

50,- = -sO-''C/-5"'=C/-'''=-|(l + ^s°°)y^ + |s^''=T/'= 

j^pklj^klj _ l^^kl^jkl _ pklyklj^ (3g) 

(37) 

Although they are unnecessary for a consistent 0(3) ex- 
pansion, the 0(3) terms for g^j and the 0(2) terms for gjk 
are displayed explicitly because they are useful for part 
of the analysis to follow. The 0(4) symbol in the ex- 
pression for (7oo serves as a reminder of the terms missing 
for a complete expansion at 0(4). Note that the met- 
ric potentials for general relativity in the chosen gauge 
arc recovered upon setting all coefficients for Lorentz vi- 
olation to zero, as expected. Note also that a nonzero 
u would merely act to scale the potentials in the above 
equations by an unobservable factor (1 + m). 

The properties of this metric under spacetime trans- 
formations are induced from those of the SME action. 
As described in Sec. II A, two different kinds of space- 
time transformation can be considered: observer trans- 
formations and particle transformations. The SME is 
invariant under observer transformations, while the co- 
efficients for Lorentz violation determine both the par- 
ticle local Lorentz violation and the particle diffeomor- 
phism violation in the theory. Since the SME includes all 
observer-invariant sources of Lorentz violation, the 0(3) 
post-newtonian metric of the minimal SME given in Eqs. 
(35)- (37) must have the same observer symmetries as the 
0(3) post-newtonian metric of general relativity. 

One relevant set of transformations under which the 
metric of general relativity is covariant are the post- 
galilean transformations [33]. These generalize the 
galilean transformations under which newtonian gravity 
is covariant. They correspond to Lorentz transformations 
in the asyniptotic;ally Minkowski regime. A post-galilcan 
transformation can be regarded as the post-newtonian 
product of a global Lorentz transformation and a possi- 
ble gauge transformation applied to preserve the chosen 
post-newtonian gauge. Explicit calculation verifies that 
the 0(3) metric of the minimal SME is unchanged by an 
observer global Lorentz transformation, up to an over- 
all gauge transformation and possible effects from 0(4). 
This suggests that the post-newtonian metric of the min- 
imal SME indeed takes the same form (35)-(37) in all ob- 
server frames related by post-galilean transformations, as 
expected. 

In contrast, covariance of the minimal-SME metric 
(35)-(37) fails under a particle post-galilean transforma- 
tion, despite the freedom to perform gauge transforma- 
tions. This behavior can be traced to the invariance (as 
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opposed to covariancc) of the coefficients for Lorentz vio- 
lation appearing in Eqs. (35)-(37) under a particle post- 
galilean transformation, which is a standard feature of 
vacuum values arising from spontaneous Lorentz viola- 
tion. The metric of the minimal SME therefore breaks 
the particle post-galilean symmetry of ordinary general 
relativity. 



B. Dynamics 

Many analyses of experimental tests involve the equa- 
tions of motion of the gravitating sources. In particu- 
lar, the many-body equations of motion for a system of 
massive bodies in the presence of Lorentz violation are 
necessary for the tests wc consider in this work. Here, 
we outline the description of massive bodies as perfect 
fluids and obtain the equations of motion and action for 
the many-body dynamics in the presence of Lorentz vio- 
lation. 



1. The post-newtonian perfect fluid 

Consider first the description of each massive body. 

Adopting standard assumptions [30] , wc suppose the ba- 
sic properties of each body arc adequately described by 
the usual energy-momentum tensor [Tm)^^ for a perfect 
fluid. Given the fluid element four- velocity m^, the mass 
density p, the internal energy 11, and the pressure p, the 
energy-momentum tensor is 

{TmY'' = {p + p\l + p)u''u'' + pg^^''. (38) 

The four equations of motion for the perfect fluid are 

D^iTMT'' = 0. (39) 

Note that the construction of the linearized effective Ein- 
stciin equations (23) in Sec. II B ensures that this equation 
is satisfied in our context. 

To proceed, we separate the temporal and spatial com- 
ponents of Eq. (39) and expand the results in a post- 
newtonian series using the metric of the minimal SME 
given in Eqs. (35)-(37), together with the associated 
Christoffel symbols. As usual, it is convenient to de- 
fine a special fluid density p* that satisfies the continuity 
equation 

dop* +d,{p*v^) = 0. (40) 
Explicitly, we have 

P* = p{-9Y'^u\ (41) 
where g is the determinant of the post-newtonian metric. 



For the temporal component = 0, we find 

= do[p*{i + n + \v^vi + u)] 

+dj [p* (1 + n + \v^v^ + uy + pv^] 

-p*doU - 2p*v^djU 

+ i9o[p*(3s"°C/ + s^''= [/■'■'=)] 

+ i9,[p*(3s°"C/ + s"C/'=>^] 

-^p*do{3-f°U + s^''W'') 

-p*v^dj{3-f'^U + s''^U''^). (42) 

The first four terms of this equation reproduce the usual 
generalized Eulcr equations as expressed in the post- 
newtonian approximation. The terms involving 's^'^ rep- 
resent the leading corrections due to Lorentz violation. 

The key effects on the perfect fluid due to Lorentz vio- 
lation arise in the spatial components v = j oi Eq. (39). 
These three equations can be rewritten in terms of p* 
and then simplified by using the continuity equation (40). 
The result is an expression describing the acceleration of 
the fluid elements. It can be written in the form 

P*^+4 + ^iy=0, (43) 

where A'^ is the usual set of terms arising in general 
relativity and A-'j^y contains terms that violate Lorentz 
symmetry. 

For completeness, we keep terms in A'-^ to 0(4) and 
terms in A-^y to 0(3) in the post-newtonian expansion. 
This choice preserves the dominant corrections to the 
perfect fluid equations of motion due to Lorentz viola- 
tion. Explicitly, we find that is given by 

= -p*d,U + dj[p{l + W)] 
-(5jp)(n+i^;V+p/p*) 

-vi{p*doU -dop) 
+Ap*d{Uv^ -V^)/dt 
+\p*do{Vi -W^) 
+Ap*v''djV'' - p*djf^ 
-p^idjUXv^v^ + Sp/p*), (44) 

where $ is a perfect-fluid potential at 0(4) given by 

+ n{^,t)+2,p{x',t)/p{x',t)]. (45) 
Similarly, the Lorentz-violating piece is found to be 

-p*dois'°U + s°''W'') 
-p*v''dkis'°U + s°^W^) 
+p*v^dj{-^*'U + ^^U*'^) 
+2p*djS°^V^. (46) 
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We note in passing that Eqs. (40) and (43) can be used 
to demonstrate that the standard formula 



dn 



= -Pdkv'' 



(47) 



remains valid in the presence of Lorentz violation. 



2. Many-body dynamics 

As a specific relevant application of these results, we 
seek the equations of motion for a system of massive bod- 
ies described as a perfect fiuid. Adopting the general 
techniques of Ref . [30] , the perfect fluid can be separated 
into distinct self-gravitating clumps of matter, and the 
equations (43) can then be appropriately integrated to 
yield the coordinate acceleration of each body. 

The explicit calculation of the acceleration requires the 
introduction of various kincmatical quantities for each 
body. For our purposes, it suffices to define the conserved 
mass of the ath body as 



ma = j d^x'p*{t,x'), 



(48) 



where p* is specified by Eq. (41). Applying the continuity 

properties of p*, we can introduce the position x^, the 
velocity vl, and the acceleration of the ath body as 

4 = m-' J d^xp*{t,x)x^, (49) 

vi = ^=m-'Jd'xp*{t,x)^, (50) 



fxi 
~dfi 



'^=m-' J d^xp*it,x) 



df^ ■ 



(51) 



The task is then to insert Eq. (43) into Eq. (51) and 
integrate. 

To perform the integration, the metric potentials are 
split into separate contributions from each body using 
the definitions (49) and (50). We also expand each po- 
tential in a multipole series. Some of the resulting terms 
involve tidal contributions from the finite size of each 
body, while others involve integrals over each individual 
body. In conventional general relativity, the latter van- 
ish if equilibrium conditions are imposed for each body. 
In the context of general metric expansions it is known 
that some parameters may introduce self accelerations 
[30], which would represent a violation of the gravita- 
tional weak equivalence principle. However, for the grav- 
itational sector of the minimal SME, we find that no 
self contributions arise to the acceleration at post- 
newtonian order 0(3) once standard equilibrium condi- 
tions for each body are imposed. 

For many applications in subsequent sections, it suf- 
fices to disregard both the tidal forces across each body 
and any higher multipole moments. This corresponds to 



taking a point-particle limit. In this limit, some calcu- 
lation reveals that the coordinate acceleration of the ath 
body due to N other bodies is given by 



fxi 
df^ 



= -(i + i^°°)E%^ 



bjta 



' ab 



' ab 



b^a 



.-J)k Sr^ ^''^bVh'''ab'''ab'''ab 



+35°'= E 
b^a 

-I-.... 



' ab 



(52) 



The ellipses represent additional acceleration terms aris- 
ing at 0(4) in the post-newtonian expansion. These 
higher-order terms include both corrections from general 
relativity arising via Eq. (44) and 0(4) effects controlled 
by SME coefficients for Lorentz violation. 

Note that these point-mass equations can be general- 
ized to include tidal forces and multipole moments using 
standard techniques. As an additional check, we verified 
the above result is also obtained by assuming that each 
body travels along a geodesic determined by the pres- 
ence of all other bodies. The geodesic equation for the 
ath body then yields Eq. (52). 

The coordinate acceleration (52) can also be derived as 
the equation of motion from an effective nonrelativistic 
action principle for the N bodies. 



dtL. 



(53) 



A calculation reveals that the associated lagrangian L 
can be written explicitly as 



ah 



Tab 



-IT. 



ah 



Tab 



(54) 



where it is understood that the summations omit b = a. 
Note that the 0(2) potential has a form similar to that 
arising in Lorentz- violating electrostatics [34, 35]. 

Using standard techniques, one can show that there 
are conserved energy and momenta for this system of 
point masses. This follows from the temporal and spa- 
tial translational invariance of the lagrangian L, which 
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in turn is a consequence of the choice d^s"^ = in Eq. 
(11). 

As an illustrative example, consider the simple case of 
two point masses m\ and in the limit where only 0(3) 
terms are considered. The conserved hamiltonian can be 
written 



H 



Pi 



p2 



Gmim2 



2mi 2m2 7 
3Gs°-'(m2P/ +miP|) 



(1 



+ 



+ 



2r 



(55) 



In this expression, r = fi — r2 is the relative separation 
of the two masses, and the canonical momenta are given 
by Pn = dL/drn for each mass. It turns out that the 
total conserved momentum is the sum of the individual 
canonical momenta. Explicitly, we find 



3Gmim2S 



■oj 



Gm\m2S^^r^r^ 



, (56) 



where the first term is the usual newtonian center-of-mass 
momentum P/v — mivi + 'm2V2- Defining M = mi +m2, 
we can adopt V = Pn /M as the net velocity of the bound 
system. 

Further insight can bo gained by considering the time 
average of the total conserved momentum for the case 
where the two masses are executing periodic bound mo- 
tion. If we keep only results at leading order in coeffi- 
cients for Lorentz violation, the trajectory r for an ellip- 
tical orbit can be used in Eq. (56). Averaging over one 
orbit then gives the mean condition 



where 



Gmim2 



< P >= M <V > +a, 



(57) 



[Se's""^ +e{e)spP^ + {l-ey^'e{e)sQQ^)] 

(58) 

is a constant vector. In the latter equation, a is the 
semimajor axis of the ellipse, e{e) is a function of its ec- 
centricity e, P points towards the periastron, while Q is 
a perpendicular vector determining the plane of the or- 
bit. The explicit forms of £(e), P, and Q are irrelevant 
for present purposes, but the interested reader can find 
them in Eqs. (166) and (169) of Sec. VE. Here, we re- 
mark that Eq. (57) has a parallel in the fermion sector 
of the minimal SME. For a single fermion with only a 
nonzero coefficient for Lorentz violation of the type, 
the nonrelativistic limit of the motion yields a momentum 
having the same form as that in Eq. (57) (see Eq. (32) of 
the first of Refs. [6]). Note also that the conservation of 
P and the constancy of a imply that the measured mean 
net velocity <V>o{ the two-body system is constant, a 
result consistent with the gravitational weak equivalence 
principle. 



C. Connection to other post-newtonian metrics 

This subsection examines the relationship between the 
metric of the pure-gravity sector of the minimal SME, 
obtained under assumptions (i)-(v) of Sec. II B, and some 
existing post-newtonian metrics. We focus here on two 
popular cases, the parametrized post-newtonian (PPN) 
formalism [29] and the anisotropic-universe model [36]. 

The philosophies of the SME and these two cases are 
different. The SME begins with an observer-invariant 
action constructed to incorporate all known physics and 
fields through the Standard Model and general relativity. 
It categorizes particle local Lorentz violations according 
to operator structure. Within this approach, the dom- 
inant Lorentz-violating effects in the pure-gravity sec- 
tor are controlled by 20 independent components in the 
traceless coefficients u, s^^", J'^'^^" _ Nq assumptions about 
post-newtonian physics are required a priori. In con- 
trast, the PPN formalism is based directly on the expan- 
sion of the metric in a post-newtonian series. It assumes 
isotropy in a special frame, and the primary terms are 
chosen based on simplicity of the source potentials. The 
corresponding effects involve 10 PPN parameters. The 
anisotropic-universe model is again different: it develops 
an effective TV-body classical point-particle lagrangian, 
with leading-order effects controlled by 11 parameters. 
These differences in philosophy and methodology mean 
that a complete match cannot be expected between the 
corresponding post-newtonian metrics. However, partial 
matches exist, as discussed below. 



1. PPN formalism 

Since the 0(4) terms in goo in the minimal-SME metric 
(35) remain unknown, while many of the 10 PPN param- 
eters appear only at 0(4) in goo, a detailed comparison 
of the pure-gravity sector of the minimal SME and the 
PPN formalism is infeasible at present. However, the ba- 
sic relationship can be extracted via careful matching of 
the known metrics in different frames, as we demonstrate 
next. 

We first consider the situation in an observer frame 

that is comoving with the expansion of the universe, often 
called the universe rest frame. In this frame, the PPN 
metric has a special form. Adopting the gauge (34) and 
keeping terms at 0(3), the PPN metric is found to be 



goo = -l + 2f7 + 0(4), 
90j = -i(67+l + i«i)X^^-i(3 
gjk = P>'[l + {3^-l)U]-{j-l)U- 



27+iai)VF^ 



jk 



(59) 



Of the 10 parameters in the PPN formalism, only two 
appear at 0(3) in this gauge and this frame. The reader 
is cautioned that the commonly used form of the PPN 
metric in the standard PPN gauge differs from that in 
Eq. (59) by virtue of the gauge choice (34). 
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The PPN formalism assumes that physics is isotropic 
in the universe rest frame. In contrast, the SME aUows 
for anisotropies in this frame. To compare the corre- 
sponding post-newtonian metrics in this frame therefore 
requires restricting the SME to an isotropic limit. In fact, 
the combination of isotropy and assumption (ii) of Sec. 
II B, which restricts attention to effects that are linear in 
the SME coefficients, imposes a severe restriction: of the 
19 independent SME coefficients contained in s^" and 
jkX^iv ^ only one combination is observer invariant under 
spatial rotations and hence isotropic. Within this as- 
sumption, to have any hope of matching to the PPN for- 
malism, the SME coefficients must therefore be restricted 
to the isotropic limiting form 

/sOO \ 

is°° 
^ ~ 

\ is°° / 

r^"^ = 0. (60) 

As discussed in Sec. Ill A, the coefficient u is unobserv- 
able in the present context and can be set to zero. In the 
special limit (60), the minimal-SME metric in Eqs. (35), 
(36), and (37) reduces to 

goo = -l + (2 + fs''0)[/ + O(4), 

go, = -i(7 + sO°)F■'"-i(l + |sO")W^^ 

g^j, = ^Jfe[i + (2-|sO°)C/] + |sO°[/^'^ (61) 

With these restrictions and in the universe rest frame, 
a match between the minimal-SME and the PPN met- 
rics becomes possible. The gravitational constant in the 
restricted minimal-SME metric (61) must be redefined as 

Gnew = G(l + |s°°). (62) 

The match is then 

7 = 1- |sOO. (63) 

From this match, we can conclude that the pure-gravity 
sector of the minimal SME describes many effects that 
lie outside the PPN formalism, since 18 SME coefficients 
cannot be matched in this frame. Moreover, the con- 
verse is also true. In principle, a similar match in the 
universe rest frame could be made at 0(4), where the 
minimal-SME metric in the isotropic limit can still be 
written in terms of just one coefficient s°°, while the 
PPN formalism requires 10 parameters. It follows that 
the PPN formalism, in turn describes many effects that 
lie outside the pure-gravity sector of the minimal SME, 
since 9 PPN parameters cannot be matched in this frame. 
The mismatch between the two is a consequence of the 
differing philosophies: effects that dominate at the level 
of a pure-gravity realistic action (gravitational sector of 
the minimal SME) evidently differ from those selected 



by requirements of simplicity at the level of the post- 
newtonian metric (PPN). 

Further insight about the relationship between the 
gravitational sector of the minimal SME and the PPN 
formalism can be gained by transforming to a Sun- 
centered frame comoving with the solar system. This 
frame is of direct relevance for experimental tests. More 
important in the present context, however, is that the 
transformation between the universe rest frame and 
the Sun-centered frame mixes terms at different post- 
newtonian order, which yields additional matching infor- 
mation. 

Suppose the Sun-centered frame is moving with a ve- 
locity w of magnitude \w\ \ relative to the universe 
rest frame. Conversion of a metric from one frame to the 
other can be accomplished with an observer post-galilean 
transformation. A complete transformation would re- 
quire using the transformation law of the metric, express- 
ing the potentials in the new coordinates, and transform- 
ing also the SME coefficients for Lorentz violation or the 
PPN parameters, all to the appropriate post-newtonian 
order. For the minimal-SME metric (61) in the isotropic 
limit, this procedure would include transforming s*"^, in- 
cluding also the change in the new gravitational coupling 
Gnew via its dependence on However, for some pur- 
poses it is convenient to perform only the first two of 
these steps. Indeed, in the context of the PPN formalism, 
this two-step procedure represents the standard choice 
adopted in the literature [30]. In effect, this means the 
PPN parameters appearing in the PPN metric for the 
Sun-centered frame remain expressed in the universe rest 
frame. For comparative purposes, we therefore adopt in 
this subsection a similar procedure for the isotropic limit 
of the minimal-SME metric (61). 

Explicitly, we find that the PPN metric in the Sun- 
centered frame and in the gauge (34) is given by 

500 = -I + '2.U - [ai - a2 ~ a^)w'^U 
-a^w^w'^U^^ + (2a3 - ai)w^V^ 
+{2a2 - \ax)w^{V' - W) + 0(4), 

-i(67 + 1 + iai)y^' - i(3 - 27 + ^ai)W^, 
gjk = + (37 - l)U] - (7 - l)C/^'=. (64) 

This expression depends on four of the 10 PPN parame- 
ters. It includes all terms at 0(3) and also all terms in 
500 dependent on w. Since w is 0(1), some of the latter 
are at 0(4), including those involving the two parameters 
a2 and a.^ that arc absent from the PPN metric (59) in 
the universe rest frame. The dependence on w is a key 
feature that permits a partial comparison of these 0(4) 
terms with the minimal SME and hence provides match- 
ing information that is unavailable in the universe rest 
frame. 

For the minimal SME, the isotropic limit of the pure- 
gravity sector in the Sun-centered frame and the gauge 
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(34) is found to be 

+ lw^w''s°°W'' + f s°W^' + 0(4), 

-i(7 + sO")T/-''-i(l + |s«")VK-'', 
gjk = S^'il + (2 - p'°)U] + |s°OC/^"^ (65) 

where the cocfScicnt s"*^ remains defined in tlie universe 
rest frame. Since the general results (35)-(37) for the 
minimal-SME metric take the same form in any post- 
galilean observer frame, Eq. (65) can be derived by trans- 
forming the coefficient s^" from the universe rest frame 
to the corresponding coefficients s^"^ in the Sun-centered 
frame and then substituting the results into Eqs. (35)- 
(37). Like the PPN metric (64), the expression (65) con- 
tains all terms at 0(3) along with some explicit 0(4) 
terms that depend on w. 

By rescaling the gravitational constant as in Eq. (62) 
and comparing the two post-newtonian metrics (64) and 
(65), we recover the previous matching results (63) and 
obtain two additional relationships: 



as = 0, 

7 = 1- (66) 

The vanishing of as is unsurprising. This parameter is 
always zero in semiconservative theories [30, 37], while 
assumption (i) of Sec. II B and Eq. (11) imply a constant 
asymptotic value for which in turn ensures global 
energy-momentum conservation. A more interesting is- 
sue is the generality of the condition 

ai = 4a2 (67) 

implied by Eq. (66). It turns out that this condition 
depends on assumption (iv) of Sec. II B, which imposes 
the vanishing of the independent energy-momentum con- 
tribution E'"'. The relationship between the conditions 
SA"' ^ and ai ^ 4a2 is considered further in Sec. IV 
below. 

We emphasize that all quantities in Eq. (66) are de- 
fined in the universe rest frame. For experimental tests of 
the SME, however, it is conventional to report measure- 
ments of the coefficients for Lorentz violation in the Sun- 
centered frame. The conversion between the two takes a 
simple form for the special isotropic limit involved here. 
It can be shown that the minimal-SME coefficients Sg" in 
the Sun-centered frame and the isotropic coefficient s^j 
in the universe rest frame are related by 

-Oj _ _4 j-00 

= ^{S^'' +4w^w'')s^°. (68) 




FIG. 1: Schematic relationship between the pure-gravity sec- 
tor of the minimal SME and the PPN formaUsm. 



These equations can be used to relate the results in this 
subsection to ones expressed in terms of the coefficients 
s'*'' in the Sun-centered frame. 

The relationship between the pure-gravity sector of the 
minimal SME and the PPN formalism can be represented 
as the Venn diagram in Fig. 1. The overlap region, which 
corresponds to the isotropic limit, is a one-parameter re- 
gion in the universe rest frame. This overlap region en- 
compasses only a small portion of the effects governed by 
the pure-gravity sector of the minimal SME and by the 
PPN formalism. Much experimental work has been done 
to explore the PPN parameters. However, the figure illus- 
trates that a large portion of coefficient space associated 
with dominant effects in a realistic action (SME) remains 
open for experimental exploration. We initiate the theo- 
retical investigation of the various possible experimental 
searches for these effects in Sec. V. 

We note in passing that the above matching consid- 
erations are derived for the pure-gravity sector of the 
minimal SME. However, even in the isotropic limit, the 
matter sector of the minimal SME contains numerous ad- 
ditional coefficients for Lorentz violation. Attempting a 
match between the isotropic limit of the minimal SME 
with matter and the PPN formalism would also be of 
interest but lies beyond our present scope. 



2. Anisotropic-universe model 

An approach adopting a somewhat different philosophy 
is the anisotropic-universe model [36] . This model is con- 
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servative and is based on the formulation of an effective 
classical point-particle lagrangian. Possible anisotropies 
in a preferred frame are parametrized via a set of spatial 
two-tensors and spatial vectors, and the post-newtonian 
metric is constructed. The model includes velocity- 
dependent terms and three-body interactions up to 0(4). 
The connection to the PPN formalism is discussed in Ref. 
[36]. 

A detailed match between the anisotropic-universe 
model and the classical point-particle limit of the pure- 
gravity sector of the minimal SME is impractical at 
present, since the 0(4) terms for the latter are imde- 
termined. However, some suggestive features of the rela- 
tionship can be obtained. For this purpose, it suffices to 
consider the restriction of the anisotropic-universe model 
to two-body terms at 0(3). In this limit, the lagrangian 
is 



Lav 



2 2^rnaV^ + ^}_^ (1 



ah 



Tab 



+ E ^^i^'vi + T^ri,v'^ . f1,). (69) 



Tab 



The anisotropic properties of the model at this order are 
controlled by 11 parameters, collected into one symmetric 
spatial two-tensor Vl^^ and two three- vectors ^-'^ T-' . 

To gain insight into the relationship, we can compare 
Eq. (69) with the 0(3) point-particle lagrangian (54) ob- 
tained in Sec. TUB. Consider the match in the preferred 
frame of the anisotropic-universe model. Comparison of 
Eqs. (69) and (54) gives the correspondence 



4* 12* 
2* ' 



1 -OO^jTs 



(70) 



The five parameters il^^ are determined by the SME co- 
efficients and s°°, while the two three- vectors , 
are determined by s°-^ . 

The SME is a complete effective field theory and so 
contains effects beyond any point-particle description, 
including that of the anisotropic-universe model. Even 
in the point-particle limit, it is plausible that the pure- 
gravity sector of the minimal SME describes effects out- 
side the anisotropic-universe model because the latter is 
based on two-tensors and vectors while the SME contains 
the four-tensor J'^^^^'^ ^ Nonetheless, in this limit the con- 
verse is also plausible: the anisotropic-universe model 
is likely to contain effects outside the pure-gravity sec- 
tor of the minimal SME. The point is that the match 
(70) implies the 11 parameters fi^'^, <!>■', are fixed 
by only 9 independent SME coefficients s^"^. This sug- 
gests that two extra degrees of freedom appear in the 
anisotropic-imiverse model already at 0(3). Some cau- 
tion with this interpretation may be advisable, as the 
condition = 3T^ implied by Eq. (70) in the context 
of the minimal SME arises from the requirement that 
the form of the lagrangian be observer invariant under 



post-galilean transformations. To our knowledge the ob- 
server transformation properties of the parameters in the 
anisotropic-universe model remain unexplored in the lit- 
erature, and adding this requirement may remove these 
two extra degrees of freedom. At 0(4), however, the 
full anisotropic-universe model contains additional two- 
tensors for which there are no matching additional coef- 
ficients in the pure-gravity sector of the minimal SME. 
It therefore appears likely that the correspondence be- 
tween the two approaches is again one of partial overlap. 
It is conceivable that effects from the SME matter sector 
or from nonminimal SME terms could produce a more 
complete correspondence. 



IV. ILLUSTRATION: BUMBLEBEE MODELS 

The analysis presented in Sees. II and III applies to 
all theories with terms for Lorentz violation that can 
be matched to the general form (5). An exploration of 
its implications for experiments is undertaken in Sec. V. 
Here, we first take a short detour to provide a practical 
illustration of the general methodology and to illuminate 
the role of the five assumptions adopted in the lineariza- 
tion procedure of Sec. II B. For this purpose, we consider 
a specific class of theories, the bumblebee models. Note, 
however, that the material in this section is inessential 
for the subsequent analysis of experiments, which is in- 
dependent of specific models. The reader can therefore 
proceed directly to Sec. V at this stage if desired. 

Bumblebee models are theories involving a vector field 
that have spontaneous Lorentz violation induced by 
a potential V{B^). The action relevant for our present 
purposes can be written as 



/ 



1 

2k 



iaeS^-B^, + l(3e{D^B,){D^^B") 



eV{B^') + CM , 



(71) 



where a and (3 are real. In Minkowski spacetime and with 
vanishing potential V, a nonzero value of /3 introduces 
a Stiickelberg ghost. However, the ghost term with /3 
potentially nonzero is kept in the action here to illustrate 
some features of the assumptions made in Sec. II B. In 
Eq. (71), the potential V is taken to have the functional 
form 



ViB") = V{B''Bf,±b'^), 



(72) 



where 6^ is a real number. This potential induces a 

nonzero vacuum value = fe^ obeying 6^6^ = =f6^. 
The theory (71) is understood to be taken in the limit of 
Riemann spacetime, where the field-strength tensor can 
be written as B^j^ = df^Bi, — duB^. Bumblebee mod- 
els involving nonzero torsion in the more general context 
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of Riemann-Cartan spacetime are investigated in Refs. 
[4, 28]. 

Theories coupling gravity to a vector field with a vac- 
uum value have a substantial history in the literature. 
A vacuum value as a gauge choice for the photon was 
discussed by Nambu [38]. Models of the form (71) with- 
out a potential term but with a vacuum value for B'^ 
and nonzero values of a and /? were considered by Will 
and Nordtvedt [39] and by Hellings and Nordtvedt [40]. 
The idea of using a potential V to break Lorentz symme- 
try spontaneously and hence to enforce a nonzero vac- 
uum value for B^^ was introduced by Kostelecky and 
Samuel [41], who studied both the smooth 'quadratic' 
case V = —Ae(i?^i3^± 6^)^/2 and the limiting Lagrange- 
muhipHer case V = -Xe{Bi'B,, ± iP) for /3 = and 
in A'^ dimensions. The spontaneous Lorentz breaking is 
accompanied by qualitatively new features, including a 
Nambu- Goldstone sector of massless modes [28], the nec- 
essary breaking of U(l) gauge invariance [4], and impli- 
cations for the behavior of the matter sector [42], the 
photon [28], and the graviton [4, 43]. More general po- 
tentials have also been investigated, and some special 
cases with hypergeometric V turn out to be renormal- 
izablc in Minkowski spacetime [44] . The situation with a 
Lagrangc-multiplier potential for a timelike 6^ and both 
a and /? nonzero has been explored by Jacobson and col- 
laborators [45] . Related analyses have been performed in 
Refs. [46, 47]. 

In what follows, we consider the linearized and post- 

newtonian limits of the theory (71) for some special cases. 
The action for the theory contains a pure-gravity Einstein 
piece, a gravity-bumblebee coupling term controlled by 
^, several terms determining the bumblebee dynamics, 
and a matter Lagrange density. Comparison of Eq. (71) 
with the action (5) suggests a correspondence between 
the gravity-bumblebee coupling and the fields u and sf^" , 
with the latter being related to the traceless part of the 
product B'^B'^. It is therefore reasonable to expect the 
linearization analysis of the previous sections can be ap- 
plied, at least for actions with appropriate bumblebee 
dynamics. 



A. Cases with /3 = 

Consider first the situation without a ghost term, 
(3 = 0. Suppose for definiteness that a = 1. The modi- 
fied Einstein and bumblebee field equations for this case 
are given in Ref. [4] . To apply the general formalism de- 
veloped in the previous sections, we first relate the bum- 
blebee action (71) to the action (3) of the pure-gravity 
sector of the minimal SME, and then we consider the 
linearized version of the equations of motion. 

The match between the bumblebee and SME actions 
involves identifying u and s'"' as composite fields of the 



underlying bumblebee field and the metric, given by 

^ ^Q^Qu _ l^^B''B^g^0g''\ (73) 

This is a nonlinear relationship between the basic fields 
of a given model and the SME, a possibility noted in Sec. 
II A. 

Following the notation of Sec. II B, the bumblebee field 
Bi^ can be expanded around its vacuum value b'^ as 

Bi^ = ^ + 31^. (74) 

In an asympotically cartesian coordinate system, the vac- 
uum value is taken to obey 

9^6" = 0. (75) 

This ensures that assumption (i) of Sec. II B holds. 

The expansion of the SME fields u, st"", i'"^^'' about 
their vacuum values is given in Eq. (10). At leading order, 
the ^-dependent term in the bumblebee action (71) is 
reproduced by the Lorentz- violating piece (5) of the SME 
action by making the identifications 

u = \i{bc.B^ + \b^hf'h^p), 
s^"' = ^(^b^"" - \rj'"'b°'b^), 
S'^'' = ^{b'^B" + b''B'' + ^h'^'b'^b^ 

- ^v^-baB" - Iv'^'b^b^h^^), 

r^^" = 0, 

^ 0. (76) 

We see that assumption (ii) of Sec. II B holds if the com- 
bination £j3^b^ is small. 

The next step is to examine the field equations. The 
nonlinear nature of the match (73) and its dependence 
on the metric imply that the direct derivation of the lin- 
earized effective Einstein equations from the bumblebee 
theory (71) differs in detail from the derivation of the 
same equations in terms of .s'"^ and u presented in Sec. 
II B. For illustrative purposes and to confirm the results 
obtained via the general linearization process, we pursue 
the direct route here. 

Linearizing in both the metric fiuctuation /i^i^ and the 
bumblebee fluctuation Bi^ but keeping all powers of the 
vacuum value the linearized Einstein equations can 
be extracted from the results of Ref. [4]. Similarly, the 
linearized equations of motion for the bumblebee fluc- 
tuations can be obtained. For simplicity, we disregard 
the possibility of direct couplings between the bumblebee 
fields and the matter sector, which means that assump- 
tion (iii) of Sec. II B holds. The bumblebee equations of 
motion then take the form 

d^'B^, = 2V'b, - |6^J?;.., (77) 
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where the prime denotes the derivative with respect to 
the argument. Once the potential V is specified, Eq. 
(77) for the bumblebee fluctuations can be inverted using 
Fourier decomposition in momentum space. 

Consider, for example, the situation for a smooth po- 
tential V = -\e{Bi^B^ ± 62)2/2. For this case, Eq. (77) 
can be written as 

{'n^.uU - - 4A5^6,)B^ = 

(78) 

Converting to momentum space, the propagator for the 
B^^ field for both timelike and spacelike 6'^ is found to be 



{A\b°'bc,+p^)p"p^' 
4Ap2(6„p«)2 



(79) 



In this equation, p'* is the four momentum and p^ = 
p'^Pix- This propagator matches results from other anal- 
yses [44, 46]. 

The propagator (79) can be used to solve for the har- 
monic bumblebee fluctuation B^{p) in momentum space 
in terms of the metric. We find 



B^{p) 



-h''''ba 



pm^Ki 

2b°'pc, 



+ 



$,b^'R ^pt'R ^.p^'b'^baR 

' 2Kp2 + 8«;A6«Pa 2Kp%°'Pa 

^baR'"' ^p^b'^b^Rc.ff 



Kp'^ 



Kp%"Pa 



(80) 



This solution can be reconverted to position space and 
substituted into the linearized gravitational field equa- 
tions to generate the eflfective Einstein equations for hf^i,. 
At leading order in the coupling ^, these take the ex- 
pected general form (23) with the identifications in Eq. 
(76), except that the coefficient of appears as —3 
rather than unity. To match the normalization conven- 
tions adopted in Eq. (23), a rescaling of the typo dis- 
cussed in Sec. II B must be performed, setting u —u/',i. 
The calculation shows that assumption (v) in Sec. II B 
holds. Moreover, the independently conserved piece S^,^ 
of the energy-momentum tensor {Sstu)fj.i' contains only a 
trace term generating the rescaling of ^^^y This bumble- 
bee theory therefore provides an explicit illustration of a 
model that weakly violates assumption (iv) of Sec. II B. 

A similar analysis can be performed for the Lagrange- 
multiplier potential V = — Ae(i3^i3^ ± b"^), for both time- 
like and spacelike 6'^. We find that the linearized limits of 
these models are also correctly described by the general 
formalism developed in the previous sections, including 
the five assumptions of Sec. II B. 

It follows that the match (66) to the PPN formalism 
holds for all these bumblebee models in the context of the 
isotropic limit, for which b'^ is the only nonzero coefficient 



in the imiverse rest frame. Note that the condition (67) is 
valid both for the Lagrange-multiplier potential and for 
the smooth potential. In fact, the /? = model with zero 
potential term V = but a nonzero isotropic vacuum 
value for B^ also satisfies the condition (67) at leading 
order [30, 40]. 

We note in passing that the limit of zero coupling ^ im- 
plies the vanishing of all the coefficients and fluctuations 
in Eq. (76). In the pure gravity-bumblebee sector, any 
Lorentz-violating effects in the effective linearized the- 
ory must then ultimately be associated with the bum- 
blebee potential V. Within the linearization assump- 
tions we have made above, this result is compatible with 
that obtained in Ref. [28] for the effective action of the 
bumblebee Nambu-Goldstone fluctuations, for which the 
Einstein-Maxwell equations are recovered in the same 
limit. Further insight can be obtained by examining 
the bumblebee trace-reversed energy-momentum tensor 
{SB)tiv, obtained by varying the minimally coupled parts 
of the bumblebee action (71) with respect to the metric. 
This variation gives 

-2(6^6, - iry^.6"6a)y'. (81) 

Note that the composite nature of u and s^'^ means that 
this expression cannot be readily identified with any of 
the various pieces of the energy-momentum tensor intro- 
duced in Sec. II B. Using the bumblebee field equations, 
{Sb)ij.u can be expressed entirely in terms of V' and terms 
proportional to ^, whatever the chosen potential. For the 
smooth quadratic potential, insertion of the bumblebee 
modes B^^ obtained in Eq. (80) yields 

V = X{2b"B„ + b°'b^hc,0) = (82) 

This shows that all terms in {SB)fii' are proportional to 
^. thereby confirming that these modes contribute no 
Lorentz-violating effects to the behavior of h^i, in the 
limit of vanishing ^. Appropriate calculations for the 
analogous modes of B^ in the case of the Lagrange- 
multiplier potential gives the same result. 



B. Cases with (3^0 

As an example of a model that lies within the mini- 
mal SME but outside the mild assumptions of Sec. II B, 
consider the limit a = and /3 = 1 of the theory (71). 
The kinetic term of this theory is expected to include 
negative-energy contributions from the ghost term. 

Following the path adopted in Sec. IV A, we expand 
the bumblebee field about its vacuum value as in Eq. 
(74), impose the condition (75) of asymptotic cartesian 
constancy, make the identifications (76), and disregard 
bumblebee couplings to the matter sector. It then follows 
that the hmit a = 0, (3 = 1 also satisfies assumptions (i). 
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(ii) and (iii) of Sec. II B. In this limit, the bumblebee 
equations of motion become 

D'^D^B, = 2V'K - (83) 

For specific potentials, the linearized form of Eq. (83) 
can be inverted by Fourier decomposition and the prop- 
agator obtained. Here, we consider for definiteness the 
potential V = -Ae(-B^B^±&2)2/2. However, most of the 
results that follow also hold for the Lagrange-multiplier 
potential. 

Inverting and substituting into the linearized modified 
Einstein equations produces the effective equations for 
hfj^i,. For our purposes, it suffices to study the case ^ = 0. 
In the harmonic gauge, we find the perturbation h^i, is 
determined by the linearized effective equations 

-lv^,b''}/nh^f3 + {SM)^.u]. (84) 

The right-hand side of this equation fails to match the 
generic form (23) of the linearized equations derived in 
Sec. II B. In fact, the non-matter part can be regarded as 
an effective contribution to the independently conserved 
piece of the energy-momentum tensor (S'si„)^,y. We 
can therefore conclude that the ghost model with a = 0, 
/? = 1 strongly violates assumption (iv) of Sec. II B. 

The ghost nature of the bumblebee modes in this 
model implies that an exploration of propagating solu- 
tions of Eq. (84) can be expected to encounter problems 
with negative energies. Nonetheless, a post-newtonian 
expansion for the metric can be performed. For definite- 
ness, we take the effective coefficients for Lorentz viola- 
tion Kf^f^ to be small, and we adopt the isotropic-limit 
assumption that in the universe rest frame only the coef- 
ficient 6*^ is nonzero. A calculation then reveals that the 
post-newtonian metric at 0(3) in the Sun-centered frame 
and in the gauge (34) is given by 

goo = -1 + 2U - K{byw'^U + K{byw^W^W^ 

-2K,{b°fw^{V^ - W^) + 0(4), 

gjk = ^^■'=(l + 2f/), (85) 

where G has been appropriately rescaled. Comparison 

with the PPN metric (59) in the same gauge yields the 
following parameter values in this ghost model: 

ai = 0, 

a2 = -K{b°f, 

as = 0, 

7 = 1. (86) 

This model therefore fails to obey the condition (67). The 

violation of assumption (iv) evidently affects the general 
structure of the post-newtonian metric. 



In light of the results obtained above, we conjecture 
that quadratic ghost kinetic terms are associated with 
a nonzero value of that strongly violates assump- 
tion (iv) of Sec. II B and that violates the condition (67) 
at linear order in the isotropic limit. This conjecture 
is consistent with the results obtained above with the 
Lagrange-multiplier potential and with the smooth po- 
tential, both for the case /3 = and for the case /3 7^ 0. 
Moreover, studies in the context of the PPN formalism of 
various models with V = Q but a nonzero vacuum value 
for Bi^ also suggest that ghost terms are associated with 
violations of the condition (67). For example, this is true 
of the post-newtonian limit of the ghost model with van- 
ishing y = and = [39]. Similarly, inspection of 
the general case with V = Q but a ^ Q and (3 [30] 
shows that the condition (67) is satisfied at linear order 
in the PPN parameters when /? vanishes. A proof of the 
conjecture in the general context appears challenging to 
obtain but would be of definite interest. 



V. EXPERIMENTAL APPLICATIONS 

The remainder of this paper investigates various grav- 
itational experiments to determine signals for nonzero 
coefficients for Lorentz violation and to estimate the at- 
tainable sensitivities. The dominant effects of Lorentz 
violation in these experiments are associated with new- 
tonian 0{2) and post-newtonian 0(3) terms in 900 and 
with post-newtonian 0{2) terms in goj. They are con- 
trolled by combinations of the 9 coefficients s^'^. Effects 
from 0(4) terms in goo he beyond the scope of the present 
analysis. However, some effects involving 0(3) terms in 
goj or 0(2) terms in gjk are accessible by focusing on 
specific measurable signals. For example, experiments in- 
volving the classic time-delay effect on a photon passing 
near a massive body or the spin precession of a gyroscope 
in curved spacetime can achieve sensitivity to terms at 
these orders. 

We begin in Sec. VA with a discussion of our frame 
conventions and transformation properties, which are ap- 
plicable to many of the experimental scenarios considered 
below. The types of experimental constraints that might 
be deduced from prior experiments are also summarized. 
The remaining subsections treat distinct categories of ex- 
periments. Section V B examines measurements obtained 
from lunar and satellite laser ranging. Section V C con- 
siders terrestrial experiments involving gravimeter and 
laboratory tests. Orbiting gyroscopes provide another 
source of information, as described in Sec. VD. The im- 
plications for Lorentz violation of observations of binary- 
pulsar systems are discussed in Sec. VE. The sensitivies 
of the classic tests, in particular the perihelion shift and 
the time-delay effect, are considered in Sec. VF. More 
speculative applications of the present theory, for exam- 
ple to the properties of dark matter or to the Pioneer 
anomaly, are also of interest but their details lie beyond 
the scope of this work and will be considered elsewhere. 



18 




FIG. 2: Sun-centered celestial-equatorial frame. 

A. General Considerations 

1. Frame conventions and transformations 

The comparative analysis of signals for Lorentz viola- 
tion from different experiments is facilitated by adopt- 
ing a standard incrtial frame. The canonical reference 
frame for SME experimental studies in Minkowski space- 
time is a Sun-centered celestial-equatorial frame [34] . In 
the present context of post-ncwtonian gravity, the stan- 
dard inertial frame is chosen as an asymptotically iner- 
tial frame that is comoving with the rest frame of the 
solar system and that coincides with the canonical Sun- 
centered frame. The cartesian coordinates in the Sun- 
centered frame are denoted by 

= {T,X^) = {T,X,Y,Z) (87) 

and are labeled with capital Greek letters. By defini- 
tion, the Z axis is aligned with the; rotation axis of the 
Earth, while the X axis points along the direction from 
the center of the Earth to the Sun at the vernal equinox. 
The inclination of the Earth's orbit is denoted rj. The 
origin of the coordinate time T is understood to be the 
time when the Earth crosses the Sun-centered X axis at 
the vernal equinox. This standard coordinate system is 
depicted in Fig. 2. The corresponding coordinate basis 
vectors are denoted 

eH = (eT,e,7) = (9T,9j). (88) 

In the Sun-centered frame, as in any other inertial 
frame, the post-newtonian spacetime metric takes the 
form given in Eqs. (35)- (37). This metric, along with 
the point-mass equations of motion (52), forms the basis 
of our experimental studies to follow. The corresponding 
line element can be written in the form 

ds'^ = -[l-hTT{T,X) + 0{A)\dT^ 
+2hTj{T,X)dTdX'^ 
+[5jK + hjK{T, xy^dX^dX'^, (89) 

where Htt is taken to 0(3), hrj is taken to 0(3) and 
hjK is taken to 0{2). For the purposes of this work, it 



typically suffices to include contributions to the metric 
fluctuations from the Sun and the Earth. 

Of particular interest for later applications are various 
sets of orthonormal basis vectors that can be defined in 
the context of the line element (89). One useful set is 
appropriate for an observer at rest, dX'^/dT = 0, at a 
given point (T, X) in the Sun-centered frame. Denoting 
the four elements of this basis set by with /x = {t,j), 
we can write 

et = <5%[l + i/iTT(T,X) + 0(4)]eT, 

= S'^jiej - lhjK{T,X)eK] + 5' ^hTj{T,X)eT. 

(90) 

Direct calculation shows that this basis satisfies 

rfs^ = -el + (91) 

to post-newtonian order. 

Another useful set of basis vectors, appropriate for an 
observer in arbitrary motion, can be obtained from the 
basis set (90) by applying a local Lorentz transformation. 
Denoting this new set of vectors by e/i, we have 

= A-^(T)e.. (92) 

It is understood that all quantities on the right-hand side 
of this equation are to be evaluated along the observer's 
worldline, which is parametrized by proper time r. 

For the experimental applications in the present work, 
it suffices to expand the local Lorentz transformation in 
Eq. (92) in a post-newtonian series. This gives 

e~. = S'..v''R'''et + 6'-.{6''^ + ^v''v^)R^^ek. (93) 

The components (e^)- coincide with the observer's four- 
velocity u" in the Sun-centered frame. In this expression, 

is the coordinate velocity of the observer as measured 
in the frame (90), and R^'^ is an arbitrary rotation. The 
reader is cautioned that the coordinate velocities and 

typically differ at the post-newtonian level. The ex- 
plicit relationship can be derived from (90) and is found 
to be 

= 6/v^{l + i/iTT) + ^S/hjKv'' + 0(4). (94) 

2. Current bounds 

In most modern tests of local Lorentz symmetry with 
gravitational experiments, the data have been analyzed 
in the context of the PPN formalism. The discussion in 
Sec. Ill C shows there is a correspondence between the 
PPN formalism and the pure-gravity sector of the mini- 
mal SME in a special limit, so it is conceivable a priori 
that existing data analyses could directly yield partial 
information about SME coefficients. Moreover, a given 
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experiment might in fact have sensitivity to one or more 
SME coefficients even if the existing data analysis fails 
to identify it, so the possibility arises that new informa- 
tion can be extracted from available data by adopting 
the SME context. The specific SME coefficients that 
could be measured by the reanalysis of existing exper- 
iments depend on details of the experimental procedure. 
Nonetheless, a general argument can be given that pro- 
vides a crude estimate of the potential sensitivities and 
the SME coefficients that might be constrained. 

The standard experimental analysis for Lorentz viola- 
tion in gravitational experiments is based on the assumed 
existence of a preferred-frame vector w. Usually, this is 
taken to have a definite magnitude and orientation for 
the solar system and is identified with the velocity of the 
solar system with respect to the rest frame of the cosmic 
microwave background radiation [30]. The strength of 
the coupling of the vector to gravity is then determined 
by the two PPN parameters ai and a2- 

Performing a data analysis of this type under the as- 
sumption that the vector w is the source of isotropy vi- 
olations is equivalent in the SME context to supposing 
that only parallel projections of the coefficients s^"^ along 
the unit vector w contribute to any signals. This holds 
regardless of the choice of w. For the coefficients s^'^, 
there are two possible parallel projections, given by 



A more interesting constraint can be obtained by con- 
sidering implications of the observed close alignment of 
the Sun's spin axis with the angular-momentum axis of 
the planetary orbits. The key point is that the metric 
term s^'^U^'' violates angular- momentum conservation. 
As a result, a spinning massive body experiences a self- 
torque controlled by the coefficients s^'^ . For the Sun, this 
effect induces a precession of the spin axis, which over the 
lifetime of the solar system would produce a misalign- 
ment of the spin axis relative to the ecliptic plane. This 
idea was introduced by Nordtvedt and used to bound 
the PPN parameter 02 [50]. In the context of the pure- 
gravity sector of the minimal SME, the analysis is similar 
and so is not presented here. It turns out that the final 
result involves a particular combination of SME coeffi- 
cients in the Sun-centered frame given by 



SSSP 



(97) 



In this expression, 5*^ are the Sun-frame components of 

the unit vector pointing in the direction of the present 
solar spin. The coefficients are the traceless compo- 
nents of s"^^, given by 



■^JK -JK 



1 xJK-TT 



(98) 



(95) 



Under the assumption that the coefficients s^^ are small 
enough for perturbation theory to be valid, we find an 
approximate bound of sssp % 10~^^. 



Also, the coefficients f'*^^'' play no role in the analysis 
at this order, in accordance with the discussion in Sec. 
II B 4. It therefore follows that sensitivity to at most 2 of 
the 19 SME coefficients s'^" and t'^'^^'^ could be extracted 
from these types of analyses of experimental data. To- 
gether with the 10 coefficients i'^^^" ^ the 7 perpendicular 
projections given by 



s , = s —vrws 



(96) 



remain unexplored in experimental analyses to date. 

The two independent measurements of parallel projec- 
tions of SME coefficients could be obtained, for example, 
from a reanalysis of existing data from lunar laser rang- 
ing [48, 49]. A crude estimate of the sensitivity to be 
expected can be obtained by assuming that the bounds 
attained on a particular term of the PPN metric in the 
Sun-centered frame roughly parallel those that might be 
achieved for the corresponding term in the pure-gravity 
minimal-SME metric. This procedure gives s|j-' < 10~^ 

and sjj'^ < 10"'^-'^. However, these estimates are untrust- 
worthy because the substantial differences between the 
two metrics imply unknown effects on the data analysis, 
so these constraints are no more than a guide to what 
might be achieved. 



B. Luncir and Satellite Ranging 

Lunar laser ranging is among the most sensitive tests 
of gravitational physics within the solar system to date. 
Dominant orbital perturbations in the context of the 
PPN metric were obtained in Ref. [51], while oscilla- 
tions arising from a subset of the parameters for the 
anisotropic-universe model were calculated in Ref. [52]. 
Here, we obtain the dominant effects on the motion of a 
satellite orbiting the Earth that arise from nonzero coef- 
ficients s^"^ for Lorentz violation. Our results are appli- 
cable to the Earth-Moon system as well as to artificial 
satellites. 

The analysis is performed in the Sun-centered frame. 

The satellite motion is affected by the Earth, the Sim, 
and other perturbing bodies. For the analysis, it is useful 
to introduce quantities to characterize the masses, posi- 
tions, sidereal frequencies, etc., relevant for the motion. 
These are listed in Table 1. The various position vectors 
are depicted in Fig. 3. 

This subsection begins by presenting the coordinate 
acceleration of the Earth-satellite separation, which is 
the primary observable in laser-ranging experiments. To 
gain insight into the content of the resulting expressions, 
we next perform a perturbative analysis that extracts the 
dominant frequencies and corresponding amplitudes for 
oscillations driven by Lorentz violation. This analysis is 
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somewhat lengthy, and it is largely relegated to Appendix 
A, with only the primary results presented in the main 
text. Finally, we estimate the likely sensitivities attain- 
able in experimental analyses of ranging to the Moon and 
to various artificial satellites. 



FIG. 3: The Earth-satellite system in the Sun-centered frame. 



Quantity Definition 



mi 




satellite mass 


Tn2 




H/ai in mass 


Ivl - 


- mi + m2 


total Earth-satellite mass 


5m 


= 1712 ~ 'ITT'l 


Earth-satellite mass difference 


m„ 




mass of T^th perturbing body 


Mq 




Sun mass 


^1 




satellite position 


To 




mean satellite orbital distance to Earth 


4 




Earth position 


R 




mean Earth orbital distance to Sun 


R® 




Earth radius 


rJ 

' n 




position of nth perturbing body 


rJ-. 


= r( -ri = {x, y, z) 


Earth-satellite separation, with magnitude r = |fi — r2| 


Sr 




deviation of Earth-satellite distance r from the mean ro 




= {mir( + m2ri)/M 


position of Earth-satellite newtonian center of mass 






separation of newtonian center of mass and nth perturbing body 


w 




mean satellite frequency 






anomalistic satellite frequency 






nth harmonic of satellite frequency 




= ^GMq/R^ 


mean Earth orbital frequency 






mean satellite orbital velocity 


Ve 




mean Earth orbital velocity 




= v{ -vi = dr-^/dT 


relative Earth-satellite velocity 




= {miv( + m-2v()/M 


velocity of Earth-satellite newtonian center of mass 


Q = 


-- J2GMRI 


Earth quadrupole moment, with J2 « 0.001 



Table 1. Definitions for analysis of laser ranging to the Moon and to artificial satellites. 

I 



1. Earth-satellite dynamics bodies, we can write 



The coordinate acceleration a^g of the relative Earth- 
satellite separation can be obtained from Eq. (52). Al- 
lowing for perturbative effects from the Sun and other 



J _ d r'^ 

(99) 
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The first three terms in this expression are newtonian 
effects that are independent of Lorentz violation. The 
first is the newtonian acceleration for a point mass in 
the newtonian gravitational field of the Earth-satellite 
system. It is given by 

= -^r'. (100) 

In this expression and what follows, G has been rescaled 
by a factor of (1 — 3s"^"^/2), which is unobservable in 
this context. The second term is the newtonian tidal 
quadrupole term, which takes the form 

ai = E Grn.i^m^r--r^) ^^^^^ 



This expression basically represents the leading effects 
from external bodies at the newtonian level of approxi- 
mation. The third term is 



a: 



Q 



2r7 



(102) 



It represents the acceleration due to the Earth's 
quadrupole moment. 

The term a'[y in Eq. (99) contains the leading Lorentz- 
violating effects. It can be split into two pieces, 



Q^LV.ES + C^LV.tidal- 



(103) 



The first of these is given by 



a 



LV.ES — 



+ 



{SGMV^s'^^ + 2G6mv^s'^^y 



{GMV^r^s^-^ +2G6mv^r^s^-^) 

GMV^'s'^'^r^ SGMV^^s^'^r'^r'^rJ 
^3 + -5 • 

(104) 

It contains the Lorentz-violating accelerations arising 

from the Earth-satellite system alone. The term tidal 
represents the Lorentz-violating tidal accelerations due 
to the presence of other bodies. For the Earth-Moon 
case, the dominant Lorentz-violating tidal accelerations 
are due to the Sun, and we find the leading contributions 
are 



■'LV.tidal 



0|s-^"(r" -3i?^ii^r^) 



- SR^r^R 'v'^ - S^R^r^R'v'^) 



- SR^r^R^V'^ - S^R^r^R^v'^). 

(105) 



In contrast, for artificial satellite orbits the corresponding 
Lorentz-violating tidal effects from the Moon and Sun 
are suppressed by about six orders of magnitude and can 
safely be ignored. 

Finally, the ellipses in Eq. (99) represent higher-order 
tidal corrections from perturbing bodies, higher multi- 
pole corrections from the Earth's potential, and gen- 
eral relativistic corrections entering at the 0(4) post- 
newtonian level. These smaller effects are ignored in the 
analysis that follows. 



2. Perturbative expansion 

To fit laser-ranging data, Eq. (99) could be entered 
into an appropriate computer code, along with model- 
ing data for other systematics of the satellite orbit [48] . 
However, to get a sense of the types of sensitivities at- 
tainable, it is useful to study analytically the acceleration 
(99) in a perturbation scheme. Here, we focus on oscilla- 
tory corrections to the observable relative Earth-satellite 
separation r = |ri — 7^2! in the Sun-centered frame. Note 
that the boost between the Sun-centered frame and the 
Earth, where observations are in fact performed, intro- 
duces corrections to r only at 0(2) and hence corrections 
to Eq. (99) at 0(4), as can be verified by determining the 
proper distance measured by an Earth observer. Sim- 
ilarly, modifications to the deviation 6r of r from the 
mean Earth-satellite distance ro also enter only at 0(2) 
and can therefore be neglected for our analysis. See, for 
example, Ref. [51]. 

Due to its length, the explicit calculation of the 
Lorentz-violating corrections to r is relegated to Ap- 
pendix A. We focus here on the main results. In pre- 
senting them, some definitions associated with the satel- 
lite orbit are needed. These are depicted in Fig. 4. The 
orientation of the orbit is described by two angles, the 
longitude of the node a and the inclination of the orbit 
(3. It is also useful to introduce the circular orbit phase 
9, which represents the angle at time T = subtended 
between the position vector fo Sind the line of ascending 
nodes in the plane of the unperturbed orbit. Note that 
corrections to all these angles arising from the boost fac- 
tor between the Sun-centered frame and the Earth are 
negligible in the present context. 

The oscillatory radial corrections Sr arising from the 
Lorentz-violating terms in Eq. (99) take the generic form 

6r = J2[An cos(a;„T + 0„) + B„ sin(w„T + <?!.„)]. (106) 



The dominant amplitudes An and _B„ and the values of 
the corresponding phases are listed in Table 2. The 
amplitudes are taken directly from Eqs. (A38). The 
phases (/)„ are related as indicated to the circular-orbit 
phase 9 shown in Fig. 4. 
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cients for Lorentz violation as 



FIG. 4: Satellite orbital parameters in the Sun-centered 

frame. For simplicity, the Earth is shown as if it were trans- 
lated to the origin of the Sun-centered coordinates. 
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B2U, = -^s'Vo 


26" 


A2u,-^o = -wero(s" - s22)/16(w - 


wo) 26* 


B2U1-U10 = -weros^^/8(a; - wq) 


261 


= -uj{5m)vQro^'^ /M{uj - ujq) 


e 


= LLj{Sm)vorQs'^^ /M{w — wq) 





An^ = VQro{bi/b2)sQ^c 





-Bfig, = F®ro(6i/62)sne6 






Table 2. Dominant Earth-satellite range oscillations. 



-^22 
S — S 



2 r)\ -^XX 



2 o^ZZ 



sin /3 s 



8^^ = 



-02 



(cos a — sin a cos /3) s 
-|-(sin^ a — cos^ a cos^ /?) 
+2 sin a cos a(l + cos^ f3) s'^^ 
+2 sin (3 cos (3 sin a 's^^ 
— 2 sin /3 cos /3 cos a , 

— sin a cos a cos /3 
-|- cos a sm a cos p s 
-|-(cos^ a — sin^ a) cos/3s^^ 

-h cos a sin /3 s"^'^ -f sin a sin /3 , 
cos as -|- sin a s , 

— sinacos/3s -|- cos a cos /3 s 
-|-sin/3s'^^. (107) 



Different combinations of coefficients for Lorentz vio- 
lation also appear in Table 2, associated with the ampli- 
tudes of the oscillations at frequency fi®. The relevant 
combinations of these coefficients in the Sun-centered 
frame are explicitly given by 

soec = 5 (sin a sin sin /3 cos /3 

-|- sin a cos a cos rj cos^ /? 

— sin a cos r] cos a) s 

— (3 cos 11 + \ cos r] SIT? a 

+i cos a sin rj sin /3 cos (3 
-h 2 COS ?7 cos a cos /?) s 

— (3 sin 77 + i cos a cos 77 sin /3 cos /3 

+ isin2/3 sin 77)5^^, 
= -(3 I sin^ a cos^ /3 i cos^ a)s'^^ 

1 • • 2 /I T^Y 

— ^ COS a Sin a Sin /3s 
-I- 5 sin a sin /? cos /3s'^^. (108) 



The mass, position, frequency, and velocity variables 
appearing in Table 2 are defined in Table 1. The eccen- 
tricity of the orbit is denoted e, and for the analysis it 
is assiimcd to be much less than unity. The quantities 
61 and 62 arc defined in Eq. (A23) of Appendix A. The 
anomalistic frequency wq is the frequency of the natural 
eccentric oscillations. Note that wq differs from w in the 
presence of perturbing bodies and quadrupole moments. 
In the lunar case, for example, 2it/{u)o — w) » 8.9 years 
[53]. 

Some of the combinations of SME coefficients appear- 
ing in Table 2 are labeled with indices 1 and 2, which 

refer to projections on the orbital plane. The label 1 rep- 
resents projection onto the line of ascending nodes, while 
the label 2 represents projection onto the perpendicular 
direction in the orbital plane. Explicitly, the combina- 
tions are given in terms of the basic Sun-centered coeffi- 



3. Experiment 

A feature of particular potential interest for experi- 
ments is the dependence of the coefficients in Table 2 on 
the orientation of the orbital plane. To illustrate this, 
consider first an equatorial satellite, for which /3 = 0. 
The observable combinations of coefficients for Lorentz 
violation then reduce to 



S"-S22 



s = 



cos2a(s^^ - s^^) -I- 2 sin 2a s^'^, 

-isin2a(s^^-s^^) 

-Fcos2a;s'^^. (109) 



This implies that the 2a; and 2w — Wq frequency bands of 
this orbit have sensitivity to the coefficients 



■^xx -^YY 



and s"^^. If instead a polar satellite is considered, for 
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which /? = 7r/2, the same frequency bands have sensitiv- 
ity to all six coefhcients in s '^: 

Sll_s22 ^ icos2a(s^^-s^^) 

+ sin2as^^ - s^^, 
s^^ = cosas-^^ + sinas^^. (110) 

These examples show that satellites with different orbits 
can place independent bounds on coefficients for Lorentz 
violation. A similar line of reasoning shows sensitivity to 
the s"^'' coefficients can be acquired as well. 

For lunar laser ranging, the even-parity coefficients 
s"^^ can be extracted from Eqs. (107) using the appro- 
priate values of the longitude of the node and the incli- 
nation. For definitcncss and to acquire insight, we adopt 
the values a « 125° and (3 « 23.5°. However, these an- 
gles vary for the Moon due to comparatively large new- 
tonian perturbations, so some caution is needed in using 
the equations that follow. In any case, with these values 
we find the even-parity coefficients are given by 

(^11 _ s22^^ ^ 0.08(s^^ + - 2s^^) 

-0.31(s^^-s^^)-1.7s^^ 
+0.60s^^ + 0.42s^^, 
{s'\ = 0.43(3^^ - s^^) - 0.31s^^ 

-0.23s^^ - 0.33s^^. (Ill) 

Evidently, lunar laser ranging measures two independent 
combinations of the s"^^ coefficients. 

For the odd-parity coefficients s^"^, adopting the same 
values for a and /? shows that lunar laser ranging offers 
sensitivity to the combinations 





= -0.60s^^ 


+ 0.82s'^^ 


1 




= -0.53s^^ 


- 0.75s^^ 


+ OAOi^^ 


(snec)([ 


= -3.1s^^- 


l.ls^^ + 


0.094s^^, 


(sf2es)([ 


= -3.4s^^-F 


- 0.037s'^^ 


+ 0.15s'^^ 



This reveals that measuring the amplitudes at frequency 
u) and phase provides sensitivity to 2 combinations of 

the 3 independent coefficients s'^'^ . Provided one can 
also independently measure the amplitudes at frequency 
fl and zero phase, the third independent coefficient in 
s^'^ can also be accessed. 

Together the above considerations imply that lunar 
laser ranging can in principle measure at least 5 inde- 
pendent combinations of coefficients for Lorentz viola- 
tion. In particular, the odd-parity coefficients s^'^ can 
be completely determined. 

To estimate the experimental sensitivity attainable 
in lunar laser ranging, we assume ranging precision at 
the centimeter level, which has already been achieved 
[54]. With the standard lunar values ro = 3.8 x 10^ m, 
V(s = 1.0 X 10"^, w/{u> - Wo) » 120, and e = 0.055, 
we find that lunar laser ranging can attain the follow- 
ing sensitivities: parts in 10^^ on the coefficient (s^^)([ , 



parts in 10"'^'' on the combination (s^^ — *^^)(I i ^^id parts 
in 10^ on the coefficients s°^, s°^, sq^s, Sfiec- These 
sensitivities may be substantially improved with the new 
Apache Point Observatory Lunar Laser-ranging Opera- 
tion (APOLLO), which is currently under development 
in New Mexico with the goal of millimeter-ranging pre- 
cision [55]. 

Next, we consider ranging to artificial satellites. In 
practice, the realistic modeling of satellite orbits is more 
challenging than the Moon due to substantial perturba- 
tive effects. However, we show here that ranging to ar- 
tificial satellites has the potential to be particularly use- 
ful for measuring independent combinations of the s'^^ 
coefficients that are less readily accessible to lunar laser 
ranging due to the approximately fixed orientation of the 
plane of the lunar orbit. 

The relevant coefficients for a generic satellite orbit are 
given in Eqs. (107) in terms of the right ascension a and 
inclination /3. For the odd-parity coefficients, a similar 
argument applies as for the lunar case, so any one satel- 
lite can in principle measure all 3 coefficients s^'^ . Using 
reasoning similar to that for lunar laser ranging, we find 
that any one satellite can make only two independent 
measurements of combinations of the even-parity coef- 
ficients s"^^ . It follows that any 2 satellites at different 
orientations can in principle measure 3 of the 4 combina- 
tions of components of s'^^"^ to which lunar laser ranging 
is insensitive. 

The reason that only 3 of the 4 can be measured is 
that the combination s-^-^ -f -I- s^^ = s^'^ is ab- 
sent from Eqs. (107) and hence only 5 combinations of 
the 6 coefficients in s"^^ are measurable via the am- 
plitudes in Table 2 for any set of orbit orientations. 
However, higher harmonics in lu, O, and other frequen- 
cies are likely to arise from the eccentric motion of the 
Moon or satellite. Although the corresponding ampli- 
tudes would be suppressed compared to those in Table 
2, they could include terms from which the combination 
s -\-s +s =s might be measured. The numer- 
ical code used to fit the data would include these ampli- 
tudes. Note that the strength of some signals might be 
enhanced by a suitable orientation of the satellite orbit, 
since the amplitudes in Table 2 depend on the orientation 
angles a and f3 through the quantities bi and 62 defined 
in Eq. (A23). 



Sal(41it(> 


No. 


r 


o( 


m) 


J 


R<-f. 


ETALON 


2 


2.5 


X 


10^ 


65° 


[56] 


GALILEO 


30 


3.0 


X 


10^ 


56° 


[57] 


GLONASS 


3 


2.5 


X 


10^ 


65° 


[56] 


GPS 


2 


2.6 


X 


10^ 


55° 


[56] 


LAGEOS I 


1 


1.2 


X 


10^ 


110° 


[56] 


LAGEOS H 


1 


1.2 


X 


10^ 


53° 


[56] 



Table 3. Some high-orbit satellites. 



Some promising high-orbit satellite missions, both cur- 
rent and future, are listed in Table 3. Although beyond 
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our present scope, it would be interesting to determine 
whether data from ranging to any of these sateUites could 
be adapted to search for the signals in Table 2. As an 
estimate of attainable sensitivities, suppose centimeter- 
ranging sensitivity is feasible and take ro ~ 10^ m, 
- wo) « 2 X 10^, and e = 0.010. These values 
produce the following estimated sensitivities from laser 
ranging to artificial satellites: parts in 10^ on s^^ — s^^, 
parts in 10^° on s^^, parts in 10^ on s°^, s*^^, and parts 
in 10^ on sogjs, sq^^c- 

Interesting possibilities may also exist for observing 
secular changes to the orbits of ncar-Earth satellites. 
Sensitivities comparable to some of those mentioned 
above may be attainable for certain coefficients [58]. A 
detailed analysis of this situation would require incorpo- 
rating also contributions to the Earth-satellite accelera- 
tion Q!gg proportional to the spherical moment of iner- 
tia of the Earth. These contributions, which can arise 
from Lorentz- violating effects proportional to the poten- 
tial [/"'^ in Eq. (28), are small for high-orbit satellites 
or the Moon and so have been neglected in the analysis 
above, but they can be substantial for near-Earth orbits. 



C. Laboratory experiments 

This section considers some sensitive laboratory ex- 
periments on the Earth. Among ones already performed 
are gravimetcr tests, analyzed in Refs. [59, 60]. For this 
case, the basic idea is to measure the force on a test 
mass held fixed above the Earth's surface, seeking ap- 
parent variations in the locally measured value of New- 
ton's gravitational constant G. The coefficient space of 
the pure-gravity sector of the minimal SME to which this 
experiment is sensitive is explored below. Other sensitive 
existing laboratory experiments use various torsion pen- 
dula, with some already being applied to probe Lorentz 
violation in the fermion sector [11]. This section also dis- 
cusses some torsion-pendulum tests that could perform 
sensitive measurements of Lorentz symmetry in the grav- 
itational sector. 



1. Theory 

To explore these ideas, we must first determine the dy- 
namics of a test particle and establish the local accelera- 
tion at a point on the Earth's surface in the presence of 
Lorentz violation controlled by the coefficients s^"^. The 
results arc required in a suitable observer coordinate sys- 
tem for experiments. 

An observer on the surface of the Earth is accelerat- 
ing and rotating with respect to the asymptotic inertial 
space. The appropriate coordinate system is therefore 
the proper reference frame of an accelerated and rotated 
observer [61, 62]. At second order in the coordinate dis- 
tance , the metric for this coordinate system can be 



written in the form 

= -{dx^f[l + 2ah^ + {a^x^f + {Jx^f 
-oj^xh^ + Ro^^ok^'x''] 
+2dx°dx\e^^'^oj''x^ - iR'^f^yiX^'J) 
+dx^dxH6.-k - |%feA^'>)- (113) 

Here, we have introduced x*^, which coincides along the 
worldline with the observer's proper time. 

From this metric, an effective Lagrange density can 
be established for test particle dynamics in the proper 

reference frame of an accelerated and rotated observer. 
The action for a test particle of mass m takes the form 

S = -m [ djj^^ 

= j dx°L. (114) 

For present purposes, it suffices to express the lagrangian 
as a post-newtonian series. Denote the coordinate posi- 
tion of the test particle by x^ and its coordinate velocity 
by x^ , where the dot signifies derivative with respect to 
x". This velocity and the rotation velocity of the Earth's 
surface are taken to be 0(1). The lagrangian then be- 
comes 

m~^L = i(f)^ - a • f -I- ^(w X f)^ 

-hf • (w X f) - ^Ro-ok^^x'' + .... (115) 

In this expression, a is the observer acceleration, oj is 
the observer angular velocity, and the various curvature 
components are taken to contain only the dominant con- 
tributions proportional to V^f/oo- The ellipses represent 
terms that contribute at 0{4) via post-newtonian quan- 
tities and sub-dominant curvature terms. 

The dominant Lorentz-violating effects arise through 
a. To establish the local acceleration in the proper ref- 
erence frame in the presence of Lorentz violation, con- 
sider first the observer's four acceleration in an arbitrary 
post-newtonian coordinate frame. This is given by the 
standard formula 

a'' = ^+r^X«"- (116) 

We then choose the post-newtonian frame to be the Sun- 
centered frame described in Sec. V A. In this frame, the 
spacetime metric is given by Eq. (89). 

To express the local acceleration in the Sun-centered 

frame, the observer's trajectory .t" is required. This can 
be adequately described to post-newtonian order by 

a;= = e=^^'=+4. (117) 

In this equation, .xg is the worldline of the Earth and ^'^ is 
the spatial coordinate location of the observer in a frame 
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centered on the Earth, for which indices arc denoted with 
a tilde. The components of the comoving spatial basis 
in that frame are denoted e~^. They can be obtained 

from Eq. (93) with W'' = and the Earth velocity v;^ 
given by Eq. (94) with V"^ = V^. To sufHcient post- 
newtonian accuracy, the observer's coordinate location 
in the Earth frame can be taken as a vector pointing 
to the observer's location and rotating with the Earth. 
Thus, modeling the Earth as a sphere of radius R^, we 
can write 

I* = i?e (sin X cos{u^T + 4>) , sin x sin(a;eT + (^) , cos x) , 

(118) 

where x is the colatitude of the observer and ^ is a stan- 
dard phase fixing the time origin (see Rcf. [34] , Appendix 
C): (/) = W0(r0 — T), where is measured from one of 
the times when the y and Y axis coincide. 

To obtain explicitly the local acceleration, the poten- 
tials (28) must be evaluated for the Earth. They can 
be written as functions of the spatial position X in the 
Sun-centered frame, and the ones relevant here are given 
by 



U = — 



GMff 



U 



JK 



|X -fel' 
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(119) 



In these equations, Af© is the mass of the Earth. The 
quantity 1^ is the spherical moment of inertia of the 
Earth, given by the integral 



(120) 



over the volume of the Earth, where y is the vector dis- 
tance from the center of the Earth and p{y) is the density 
of the Earth. The ellipses in Eq. (119) represent neglected 
gravitomagnetic effects. 

Inserting Eq. (117) into Eq. (116) yields the follow- 
ing components of the observer's acceleration in the Sun- 
centered frame: 

= 0(3), 

-3g(iis Vg, S +i2S V(^6 j6 4 4 ) 
+5-'-.d^^' /dT^ + 0(4). (121) 



In these equations, the reference gravitational accelera- 
tion g is 

g = GMe/i?|, (122) 
and the quantities ii, i2, is are defined by 



ii = l + ^ie, 
«2 = 1 - 
^3 = 1 - i@- 



(123) 



Note that the analysis here disregards possible tidal 
changes in the Earth's shape and mass distribution 
[30, 63], which under some circumstances may enhance 
the observability of signals for Lorentz violation. 

To obtain the corresponding expressions in the ob- 
server's frame, it suffices to project Eq. (121) along the 
comoving basis vectors e/i — e^^es of the accelerated 
and rotated observer on the Earth's surface. Explicit ex- 
pressions for these basis vectors are given by Eq. (93), 
with the observer velocity vl and the rotation i?-'*' taken 
as 

vi = S/V^ + S'-.dC^/dT + 0{3), 

W'^ = R^'iiu^T). (124) 

To sufficient approximation, we may use 

Vq = V0(sinO0T, — COST? cos OqT, — sinT/cosfl^T), 

(125) 

where V© is the mean Earth orbital speed and r] is the 
inclination of the Earth's orbit, shown in Fig. 2. 

We are free to use the observer rotational invariance 
of the metric (113) to choose the observer's spatial co- 
ordinates to coincide with the standard SME conven- 
tions for an Earth laboratory (see Ref. [34], Appendix C). 
Thus, at a given point on the Earth's surface, z points 
towards the zenith, y points east, and x points south. 
With these conventions, we find the desired expression 
for the local acceleration to be 

= g{l + li^s'^^ + li2r')S^' -ghP' 
—(J^Rq (siv? x^''^ + sin x cos X'^''^) 

+giz{f'^vi + i^'vl) 

-^gihi^'Vi + i2s'^'Vi)5^' + .... (126) 

The ellipses represent newtonian tidal corrections from 
the Sun and Moon, along with 0(4) terms. The various 
projections of quantities along the local axes appearing 
in Eq. (126) are obtained using the comoving basis vec- 
tors = e" ^e-E- For example, = e-'S^^. Typically, 
these projections introduce time dependence on T. Note, 
however, that the coefficient s^-^ carries no time depen- 
dence up to 0(3) and hence represents an unobservable 
scaling. 
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2. Gravimeter tests 

By restricting attention to motion along the z axis, the 
expression (126) for the local acceleration can be used to 
obtain the dominant Lorentz-violating contributions to 
the apparent local variation 6G of Newton's gravitational 
constant, as would be measured with a gravimeter. This 
gives 

= g{l + l^is^"-^' + ^iis^^) - lu'^Rq sin^ x 

-gt4S^'vi - Sgiis^^Vi, (127) 

where 

i4 = l-3i®. (128) 



The time dependence of the apparent variations in G 
arises from the terms involving the indices z and the 
components s^"^ . It can be decomposed in frequency ac- 
cording to 

SG 

— = ^[C„ cos(w„T + (/>„) + £)„sin(w„r + ^„)]. 

n 

(129) 

Explicit calculation of the time dependence using the co- 
moving basis vectors yields the expressions for the am- 
plitudes Cn, Dn and phases (/)„ listed in Table 3. In this 
table, the angle x is the colatitude of the laboratory, and 
the other quantities are defined in the previous subsec- 
tion. 
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= \us^^ sin 2x 
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Dui = \iiS^^ sin 2x 






C2io+n = -1*4^05^^ (cos 77 - 1) sin^ x 
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D2u,+a = jIaVqs'^^ (cos rj - 1) sin^ x 




2(t) 


C2U-Q. = -\i4V(ss'^^{l + COST?) sin^ x 




24> 


D2u-n = \iiV(Qi^^{l + COST]) sin^ x 




2(j) 


Cn = Ve[s^-*^cos?7(ii4sin^x + 3ii) +' 


s^^ sin 77(3^1 + 14 cos^ x)] 





Da = -Ves^^(5i4sin2x + 3ii) 







Coj+n = li4VQs'^-^ sin2xsin?7 




<P 


D^+n = -ii4V'e[s^^(l - cos??) - s^^ 


sin T]] sin 2x 





Cu,-Q. = 1*4 V®s^^ sin 77 sin 2x (j) 
-Pgj-n = jHV® [5^^(1 + cos 77) + s^^ sin 77] sin 2x ^ 

Table 4. Gravimeter amplitudes. 



Assuming each amplitude in Table 4 can be separately 
extracted from real data, it follows that gravimeter ex- 
periments can measure up to 4 independent coefficients 
in s'^^ and all 3 coefficients in 'E^'f Moreover, gravime- 
ter experiments have attained impressive accuracies [60] 
that would translate into stringent sensitivities on the 
pure-gravity coefficients in the minimal SME. A crude es- 
timate suggests that gravimeter experiments fitting data 
to Eq. (129) could achieve sensitivity to s'^^ at parts in 
10^^ and to s^'^ at parts in 10^. Possible systcmatics that 
would require attention include tidal influences from the 
Sun and Moon. These can, for example, generate oscil- 
lations contributing to Eq. (129) at frequencies 2(jj^ and 



3. Torsion-pendulum tests 

A complementary class of laboratory tests could ex- 
ploit the anisotropy of the locally measured acceleration 
in the horizontal directions instead of the vertical one. 
Explicitly, the accelerations in the x and y directions are 
given by 

a' = -ghs'^'' - w^i?® sin x cos x 
+ghs^Wi + gi^s^^Vl, 

aP = -gias^' + giss^'vi giss^Wi (130) 

to order 0(3). In effect, the Lorentz-violating contribu- 
tions to these expressions represent slowly varying accel- 
erations in the horizontal directions as the Earth rotates 
and orbits the Sun. In contrast, the conventional cen- 
trifugal term in Eq. (130) is constant in time. 
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FIG. 5: Idealized apparatus with iV = 3 masses. 



To gain insight into the experimental possibilities for 

testing the effects contained in Eqs. (130), consider an 
idealized scenario with a disk of radius R and mass M 
suspended at its center and positioned to rotate freely 
about its center of symmetry in the x-y plane. We 
suppose that iV spherical masses, each of mass m^, 
2 = 1,2, ...,A'', are attached rigidly to the disk at ra- 
dius tq < R and initially located at angular position 0^ 
relative to the x axis. For convenience, we write the 
masses as rrij = fijui, where to is a basic mass and fj,j 
are dimensionless numbers. In what follows, we assume a 
torsion suspension and focus attention on small angular 
deviations 6 of the disk relative to its initial angular ori- 
entation. The apparatus with iV = 3 is depicted in Fig. 
5. Note that the locations of the center of mass and the 
center of suspension arc typically different. In practice, 
the technical challenge of keeping the disk level might be 
overcome in several ways, perhaps involving rigid suspen- 
sion or magnetic support, but this is a secondary issue in 
the present context and is disregarded here. 

Using Eq. (115), the effective potential energy for the 
apparatus can be constructed. Up to an overall constant, 
we find 

V = -^mroia'^CN + ay Sn + raulsin^ xC^Cn -W"^ 
~mro{a^ Sn — a^Cw + |a^© sin^ xSn2)0- (131) 

In this expression, we have introduced the following 
quantities: 



Sn — 



Sn2 = 



N 

^/Xjsin^j-, 

JV 

^/ijsin2^j. 



These quantities depend on the number N of masses 
used. Note that the effective potential (131) depends 
only on the distribution of the N masses and is indepen- 
dent of the disk mass M. 

The first part of Eq. (131) is proportional to 9^. This 
term contains contributions both from Lorentz violation 
and from the Earth's rotation. It acts as an effective 
restoring force for the disk, producing a time-dependent 
shift to the effective spring constant of the system that 
primarily affects free oscillations. The second part of 
Eq. (131) is proportional to 9, and it contains a time- 
independent contribution from the Earth's rotation along 
with a Lorcntz-violating piece. The latter produces a 
slowly varying Lorentz- violating torque t = —dV/dO on 
the disk. 

The oscillations of the system are determined by the 
second-order differential equation 



d'^9 ^ d9 



(133) 



where I is the total moment of inertia of the disk and 
masses, and 7 is the damping parameter of the torsion 
fiber. In writing this equation, we have made use of the 
relation d/dx° = {dT/dx^)d/dT = {1 + 0{2))d/dT. 

The presence of the damping term in Eq. (133) ensures 

that free oscillations vanish in the steady-state solution. 
In the absence of additional torques, the steady-state so- 
lution is given by 



sin a„ — Fn cos a„] sin(a;„T -|- /?„) 

-[En COS an + Fn sina„] cos(w„r -|- /3„)^ , 

(134) 

where 



/3„ = 2-fuJn/{ujl - ujI) 



(135) 



(132) 



and, as before, T is the time coordinate in the Sun- 
centered frame defined in Sec. VAl. Table 5 provides 
the amplitudes F„ and the phases a„ in terms of 
quantities defined above. 
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Torsion- Pendulum Amplitudes 

Amplitude Phase 

= s^^Cn sinx - ^(s^^ - s^^)Sn sin 2x 2f 

F2u> = - s^^)Cn sinx - ^s'^'^Sn sin 2x 20 

i^a, = s^^Cn cos X - s^^Sn cos 2x 

F„ = -s^^Cat cosx - s^^Sn cos2x 

E2u;+n = -5V'©s^^C]v(l -cos?7)sinx- |V©s^'*'S'Af (1 ~cos77)sin2x 20 

F2u;+n = -5V©s'^^C]v(l - COST?) sinx + jV(ss'^^Sn{1 - cosry) sin2x 20 

E2u,-n = ^V(ss^^Cn{1 + COST]) sinx + iFes^^S'Ar(l + cos??) sin2x 20 

F2u>-n = iVes^^CAr(l + COST?) sinx- |l^es^-^5'Ar(l + cosT?)sin2x 20 

Eu,+n = ^V^s^^Cat sinr/cosx + V^s^^S'Ar sinT?(i - cos^ x) 

+ |V©s^^CAr(cosT? - 1) cosx 

-FL+n = -^V(ss'^-^Cn sinT?cosx - V©s^^S'jv(l - cosT?)(i - cos^ x) (t> 

+ V(!)s'^^Sn sinr/(i - cos^ x) 

Eu-n = ^V(ss'^^Cn{1 + cost?) cosx + ^Viss^'^Cjv sin t? cosx 4> 

+ Vqs'^^Sn sinT?(i - cos^ x) 

Foj-n = -|V0S^-^CAr sin7?cosx + V(ss'^^ Sn sinT?(| - cos^ x) <P 

+ \^©s^^S'Ar(l + cosr/)(i -cos2 x) 

Eq = — i VqS-^^S'jv cos t? sin 2x + V^s^^ Sn sin t? sin 2x 

Fa = \V^s'^''Sn sin2x 0_ 



Table 5. Torsion-pendulum amplitudes. 



To estimate the attainable sensitivities to Lorentz vio- 
lation, we suppose the angular deflection of the apparatus 
can be measured at the nrad level, which in a different 
context has already been achieved with a torsion pendu- 
lum [64]. For definiteness, wc suppose that the disk ra- 
dius is ii = 10 cm and that the = 3 masses have equal 
magnitude rrij = M and arc located at radius ro « i? 
at angular locations 9i = tt/2, 62 = 47r/3, ^3 = Stt/S. 
Assuming a natural oscillation frequency luq « 10~^ Hz 
and a small damping parameter 7 « 10"'^, and taking 
is = 2/5 and m !v M, we find the apparatus could 
attain sensitivity to the coefficients s'^^ at the level of 
parts in 10^^ and to the s^'^ coefficients at parts in 10^^. 
These are idealized sensitivities, which are unlikely to be 
achieved in practice. Nonetheless, the result is of interest 
because it is about four orders of magnitude better than 
the (realistic) sensitivities that could be attained using 
current data from lunar laser ranging. 

Provided the amplitudes in Table 5 can be separately 
extracted, and assuming measiircmcnts can be made with 
two different initial orientations 9i, it follows that these 
types of torsion-pendulum experiments can measure 4 co- 
efficients in s'^^ and all 3 independent coefficients in 
s'^^ . As before, major systematics may include solar and 
lunar tidal effects. 

Modifications of the apparatus proposed above could 
lead to improved sensitivities. One theoretically simple 
possibility is to decrease the effective spring constant of 
the system. Different kinds of pendulum could also be 
used. For example, the ring of point masses around the 



disk could be replaced with a ring of holes instead. A 
more radical possibility is to change the geometry of the 
apparatus. For example, one could replace the disk and 
masses with a thin rod lying in the horizontal plane and 
hanging from a torsion fiber. If the rod is suspended off 
center by some means while able to rotate freely about 
the vertical axis of the fiber, the steady-state oscillations 
again depend on the coefficients for Lorentz violation. In 
fact, the explicit results for this scenario can be obtained 
from Eq. (134) and Table 5 by making the following re- 
placements: 

ro ^ \{L-2d), 

Cn — > COS^O) 

Sn sin6»o. (136) 

Here, L is the total length of the rod, and d is the distance 

from the end of the rod to the point of suspension. Other 
quantities still appear but arc reinterpreted: m is the 
mass of the rod, / is its moment of inertia, and 9q is 
the initial angle the rod subtends to the x axis. Again, 
the practical issue of suspending a rod off center could 
conceivably be addressed with a rigid support system or 
magnetic levitation. The former method might allow an 
increase in sensitivity by reducing the effective torsion 
fiber frequency ujq of the system. In any case, it would 
be of definite interest to investigate various experiments 
of this type. 

Another potentially important issue is whether any of 
the existing torsion-pendulum experiments have sensitiv- 
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ity to Lorentz-violating effects. Consider, for example, an 
experiment designed to search, for deviations from New- 
ton's law of gravity at sub-millimeter distances [64] . The 
basic component of this experiment is a pair of disks po- 
sitioned one above the other and with symmetric rings of 
holes in each. The upper disk is suspended as a torsion 
pendulum, in a manner analogous to that depicted in Fig. 
5, while the lower one rotates uniformly. A key effect of 
the holes is to produce a time- varying torque on the pen- 
dulum that varies with the vertical separation between 
the two disks. Measurements of the torques at the rota- 
tion frequency and its harmonics offers high sensitivity to 
short-distance deviations from Newton's law of gravity. 
The question of interest here is whether the presence of 
Lorentz violation would change the predicted newtonian 
torque in an observable way. A detailed analysis of this 
experiment is involved and lies beyond our present scope. 
However, some insight can be gained by using a simpli- 
fied model that treats the rings of holes as two rings of 
point masses and comparable radii positioned one above 
the other, with the lower ring rotating. 

Consider first the case of two point masses rni and m2 
at coordinate locations xi and X2. The modified newto- 
nian potential V predicted by the pure-gravity sector of 
the minimal SME is 

where x = {xi — X2)/|xi — X2\- The non-newtonian ef- 
fects are controlled by the SME coefficients for Lorentz 
violation s-'*. Although Eq. (137) maintains the inverse- 
distance behavior of the usual newtonian potential, the 
associated force is misaligned relative to the vector x be- 
tween the two point masses. It is therefore conceivable 
that an unconventional vertical dependence of the torque 
might arise in the experiment discussed above. To inves- 
tigate this rigorously, the result (137) would be inserted 
into the computer code that models the experiment and 
determines the predicted torque [64]. 

In the context of the simplified model involving two 
rings of point masses, the difference between the newto- 
nian and Lorentz-violating torques can be explored as- 
suming a single nonzero coefficient s-''^ in Eq. (137). We 
examined the Fourier components of the time-varying 
torque on the upper disk at harmonics of the rotational 
frequency of the bottom disk. The amplitudes of these 
Fourier components depend sensitively on the number 
of masses used, on the geometry of the experiment, and 
on the particular harmonic considered. This indicates 
that a careful analysis involving a detailed model of the 
actual experiment is necessary for definitive predictions. 
However, the results suggest that Lorentz-violating ef- 
fects from certain coefficients s^*^ of order 10^'^ to 10"^ 
could be discernable at a separation of about 100 microns 
in experiments of this type, at the currently attainable 
torque sensitivities of about lO"""^^ Nm. Other experi- 
ments studying deviations from short-range newtonian 



gravity [65, 66] may also be sensitive to these types of 
effects. 



D. Gyroscope Experiment 

In curved spacetime, the spin of a freely falling test 
body typically precesses relative to asymptotically flat 
spacetime [67, 68]. The primary relativistic effects for 
a gyroscope orbiting the Earth are the geodetic or dc 
Sitter precession about an axis perpendicular to the 
orbit and the gravitomagnetic frame-dragging or Lens- 
Thirring precession about the spin axis of the Earth. In 
the present context of the post-newtonian metric for the 
pure-gravity sector of the minimal SME given in Eqs. 
(35)- (37), we find that precession also occurs due to 
Lorentz violation. In particular, there can be precession 
of the spin axis about a direction perpendicular to both 
the spin axis of the Earth and the angular-momentum 
axis of the orbit. 



1. Theory 

To describe the motion of an orbiting gyroscope, 
we adopt standard assumptions [62] and use the 
Fermi- Walker transport equation given in coordinate- 
independent form by 

VuS = u(a-S), (138) 

where u is the four-velocity of the gyroscope, S is the 
spin four-vector and a is the four-acceleration. The dot 
denotes the inner product with the metric two-tensor g. 

To obtain the equation of motion for the spin vector S, 
we first construct a local frame that is comoving with 
S but is nonrotating with respect to the fixed stars in the 
asymptotically flat spacetime. The construction of this 
frame proceeds as in Sec. VA, using the post-newtonian 
metric in Eqs. (89). The basis can be taken from Eq. (93) 
by setting the rotation terms to zero. Using Eq. (138), 
we can then find the proper-time derivative of S • e--. It 
is given by 

dS ■ e-. 

^ = S-(Vuej.). (139) 

This represents the proper time rate of change of the 
spatial spin vector as measured in a comoving but non- 
rotating frame. 

To display the dependence of the spin precession on the 
underlying geometry, which includes Lorentz-violating 
terms in the metric, the right-hand side of Eq. (139) can 
be expressed in the Sun-centered frame. This requires 
expressing the components of S and VuS-- to the appro- 
priate post-newtonian order. Of particular use for this 
purpose are the Sun-centered frame components of e-.. 
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which are given by 

e^, = 5-'.[v'{l + hTT + W) 

+ lhjKv'' + hTj]+0{A), 

e'-. = S'^.idjK - ^hjK + k^^v''] + 0{4). (140) 
We can then estabhsh that 





= v-^ 3-^ + 0(2), 










(Vue.)^ 


= S'.a' + 0{3), 




(Vue^)^ 


= 5'.[\v''d.jhTT-\v'dKhTT + 






+ d^jhK]T + V^d[jhK]L] + 0(4) 


• (141) 



Combining Eqs. (140) and (141) via the expression 

S-(Vue3.) = 5=(Vue3.)"5Hn (142) 

yields the general expression for the rate of change of 
the spatial spin vector of a gyroscope, vaUd to post- 
newtonian 0(3). It is given by 

+ dljhK]T+v'^dljhK]L)- (143) 

In the Hmit of the usual general-relativistic metric, this 
agrees with existing results [30, 62] . 

2. Gyroscope in Earth orbit 

In this subsection, wc specialize to a gyroscope in a 
near-circular orbit around the Earth. For this case, the 
acceleration term vanishes. 

The first step is to determine the contributions to the 
metric potentials (28) from the various relevant bodies in 
the solar system. For simplicity, we focus on the poten- 
tials due to the Earth, neglecting any contributions from 
the Sun, Moon, or other sources, and we take the Earth 
as a perfect sphere rotating uniformly at frequency u) and 
with spherical moment of inertia give by Eq. (120). 
For convenience in what follows, we introduce a quantity 
% related to I® along with corresponding quantities if^a) 
defined for any real number a by the equations 



He) 



= l + ctie) 



(144) 



where ro is the mean orbital distance to the gyroscope. 

The explicit form is required for five of the six 
types of metric potentials in Eq. (28). Some general- 
relativistic contributions from the potential W'^ involv- 
ing the Earth's angular momentum = 2Iquj^ /'i can 



be converted into ones involving V'' by using the identity 
(30). Also, contributions from the potential y^^^ can 
be converted into ones involving only the other potentials 
by virtue of the antisymmetrization in Eq. (143) and the 
identity (33). The five necessary Earth potentials can 
be written as functions of the spatial position X in the 
Sun-centered frame, given by 



U 



GMq{X - RY {X - R)^ 



\X-R[ 



Gla 



3\X-R\5 



^[3{X - Ry {X - Rf 



V = 



Ge-'^^J^{X-Ry 



■^JKL 



\X-R\ 2\X-R\3 
ViU'^'^ + . . . , (145) 



where R = RR is the Earth's location in the Sun-centered 
frame and is its velocity. Explicit expressions for 
these quantities are given in Eq. (A9) of Appendix A 
and in the associated discussion. 

The ellipses in Eqs. (145) represent omitted pieces of 
the potentials arising from the purely rotational com- 
ponent of the Earth's motion, which generate gravito- 
magnetic effects involving Lorentz violation. A crude 
estimate reveals that these gravitomagnetic terms are 
suppressed relative to the associated geodetic terms dis- 
played in Eqs. (145) by about two orders of magnitude. 
As a result, although gravitomagnetic frame-dragging 
terms are crucial for tests of conventional general relativ- 
ity [30, 62], they can be neglected for studies of the dom- 
inant Lorentz- violating contributions to the spin preces- 
sion. For simplicity, the frame-dragging precession effects 
involving Lorentz violation are omitted in what follows. 

For use in Eq. (143), these potentials must be eval- 
uated along the worldline of the gyroscope. At post- 
newtonian order, this worldline is adequately described 
by the expressions for an arbitrary satellite orbit given 
in Appendix A. Thus, the value of {X — R) along the 
worldline can be written as 



{X - -R)gyro = Top, 



(146) 



where p is given by Eq. (A6). To post-newtonian order, 
the gyroscope velocity can be taken as v^^^ — ^® + 
, where fg is the mean orbital velocity defined by wg = 
d{rop)/dT = vqt. The normal to the plane of the orbit 
is a- = p X f. 

Inserting the potentials (145) into Eq. (143) and eval- 
uating along the worldline yields a result that can be 
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separated into two pieces, 



dT 



dSP dS? 
j_ J i_ 

dT ' 



dT 



(147) 



The first piece is the standard result from general rela- 
tivity, given by 



d5f 
j_ 

dT 



= 5<5^55^[(-34Vi + i4V) 



(148) 



where the reference gravitational acceleration g is now 

g = GM^/rl (149) 

The second piece is proportional to the coefficients for 
Lorentz violation s"^^, and it takes the form 



dS? 



2^^[^p-^l+4z(_i)S^[^t;o^lp^ 



+4^p^l(-P(3/5)^^^ 



5«(3/5)S'' - fZ(_5/3)S-^^P^/5^) 



■LM -M 



(150) 

Among the effects described by the result (147) arc os- 
cillations of the orientation of the gyroscope spin vector 
that occur at frequencies ui, uj± Qq, and higher harmon- 
ics. At the frequency cu, these oscillations change the 
angular orientation of a gyroscope by roughly 69 w 2TrvQ . 
For a gyroscope in low Earth orbit, this is 56 ^ 8x 10~^ 
arcseconds, well below observable levels for existing sen- 
sitivity to nonsecular changes. 

We therefore focus instead on the dominant secular 
contributions to the motion of the gyroscope spin. Aver- 
aging various key quantities over an orbital period yields 
the following results: 



2ro(d -£7 £7 ), 



< rMvir^ > « 0, 
< riV^ > « 0, 

< > = 0. (151) 

In these equations, the vector o"^ is given explicitly as 

a"^ = (sinasin/3, — cosasin/3, cos/3), (152) 

where the angles a and /3 are those defined in Sec. VB 
and Fig. 4. Note that the averaging process eliminates 
the terms involving VA- In the conventional case, such 



terms can be shown to be unobservable in principle by a 
consideration of the physical precession referenced to the 
fixed stars [30, 69]. It is an open question whether this 
stronger result remains true in the presence of Lorentz 
violation. 

Combining the above results yields the vector equation 
for the secular evolution of the gyroscope spin. It can be 
written in the form 



d5 
dt 



gvQ^ X S. 



(153) 



The secular precession frequency Q. is comprised of two 
pieces. 



(154) 



The first term fl^, contains precession due to conventional 
effects in general relativity, and it is given by 
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The second term 17^ contains contributions from the co- 
efficients for Lorentz violation s"^^, and it has the form 



= 8V*(-i/3)^ -«(-5/3)S a a )a 



+ 4«(-3/5)S cr . 



(156) 



In the limit of no Lorentz violation, Eq. (154) reduces 
as expected to the conventional result (155) of gen- 
eral relativity involving the geodetic precession arising 
from the spacetime curvature near the Earth and the 
frame-dragging precession arising from the rotation of 
the Earth. 



3. Gravity Probe B 

Gravity Probe B (GPB) [70] is an orbiting-gyroscope 
experiment that has potential sensitivity to certain com- 
binations of coefficients contained in Eq. (156). Our gen- 
eral results (153) for the spin precession can be special- 
ized to this experiment. In terms of the angles in Fig. 4, 
the GPB orbit is polar with a = 163° and /3 = 90°. One 
consequence of the polar orbit is that the normal vector 
(T"' to the orbital plane appearing in Eq. (156) reduces to 
the simpler form 



a'^ = (sina, — cosa, 0). 



(157) 



The terms in Eq. (153) proportional to axS act to pre- 
cess the gyroscope spin in the plane of the orbit. They 
include the geodetic precession in Eq. (155) predicted 
by general relativity, along with modifications due to 
Lorentz violation. The magnitude of the precession is 
determined by the projection a ■ Q.. The explicit form of 
the relevant projection of Q. arising due to Lorentz viola- 
tion is found to be 
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The terms in Eq. (153) proportional to J x S, where 
J is the Earth's angular momentum, cause precession of 
the gyroscope spin in the plane perpendicular to J. Since 
J lies along Z, this is the XY plane of Fig. 4. These 
terms include the frame-dragging effect in Eq. (155) pre- 
dicted by general relativity, together with effects due to 
Lorentz violation. These have magnitude determined by 
the projection Z-(l. For the contributions due to Lorentz 
violation, we find the explicit form 

^^f = |i(-3/5)S^'^^'^. (159) 

In addition to the above, however, Lorentz violation 
induces a third type of precession that is qualitatively 
different from those predicted by general relativity. This 
is a precession of the gyroscope spin in the plane with 
normal h along a x J. With the definitions shown in 
Fig. 4, this normal is given by n = (cos a, sin a, 0). The 
magnitude of the associated precession is given by n • fi. 
We obtain 

^ = P(-3/5)S^''«^a^. (160) 

The GPB experiment is expected to attain sensitivity 
to secular angular precessions of the order of 5 x 10^'* 
arcseconds per year. If the spin precession in the three 
orthogonal directions can be disentangled, GPB could 
produce the first measurements of three combinations of 
the s'^^ coefficients, including one involving precession 
about the h direction. The attainable sensitivity to s'^^ 
coefficients is expected to be at the 10~^ level. 

E. Binary pulsars 

Binary pulsars form a useful testing ground for general 
relativity [71, 72]. This subsection examines the possi- 
bility of probing SME coefficients for Lorentz violation 
using pulsar timing data [73]. We find that Lorentz vi- 
olation induces secular variations in the orbital elements 
of binary-pulsar systems, along with modifications to the 
standard pulsar timing formula. The discussion is sep- 
arated into five parts: one establishing the basic frame- 
work and assumptions, a second presenting the secular 
changes, a third describing the timing formula, a fourth 
determining some experimental effects, and a final part 
explaining the transformation between the frame used to 
study the properties of the binary pulsar and the stan- 
dard Sun-centered frame. 

The reader is cautioned that the notation for the 
binary-pulsar systems in this subsection is chosen to 
match the literature and differs in certain respects from 
the notation adopted elsewhere in the paper. This nota- 
tion is described in Sec. VEl. 

1. Basics 

The general problem of describing the dynamics of bi- 
nary pulsars is complicated by the strong-field gravita- 



tional properties of the spacetime regions close to the 
pulsar and its companion. Nonetheless, provided the typ- 
ical radius of the pulsar and companion are much smaller 
than the average orbital distance between them, the or- 
bital dynamics can be modeled with an effective post- 
newtonian description. 

One standard approach is to construct a modified 
Einstein-Infeld-Hoffman (EIH) lagrangian describing the 
post-newtonian dynamics of the system [30]. In princi- 
ple, an analysis of the effects of Lorentz violation could 
be performed via a modified EIH approach based on the 
action of the pure-gravity sector of the minimal SME. 
This would involve generalizing Eq. (54) to allow for 
the usual EIH effects, along with modifications to the 
Lorcntz-violating behavior arising from the structure of 
the pulsar and its companion. A detailed analysis along 
these lines would be of interest but lies beyond the scope 
of the present work. 

In fact, a simpler approach suffices to extract key fea- 
tures arising from Lorentz violation. Instead of a full 
EIH- type treatment, we adopt in what follows a sim- 
plified picture and consider only the effects arising in a 
point-mass approximation for each of the two bodies in 
the system. To retain some generality and an EIH flavor 
in the treatment, we interpret the gravitational constant, 
the masses, and the coefficients for Lorentz violation as 
effective quantities [74]. However, the results disregard 
details of possible strong-field gravitational effects or con- 
sequences of nonzero multipole moments and tidal forces 
across the pulsar and companion. The reader is cau- 
tioned that the corresponding conventional effects from 
newtonian theory can be crucial for studies of binary pul- 
sars, so they should be included in a complete analysis. 
Nonetheless, the results obtained here offer a good base- 
line for the sensitivities to Lorentz violation that might 
be achieved through observations of binary pulsars. 

The analysis that follows makes use of the standard 
keplerian characterization of an elliptical orbit [75] . Some 
of the key quantities are displayed in Fig. 6. In this figure, 
the point P is the location of the pulsar at time t, and 
it lies a coordinate distance r from the focus F of the 
elliptical orbit. The line segments CA and CB are the 
semimajor axes, each of length a. The eccentric anomaly 
E is defined as the angle ACQ, where Q is the point on 
the enclosing circle with the same abcissa as P. The true 
anomaly / is defined as the angle AFP, while w is the 
angle subtended between the line of ascending nodes FN 
and the major axis of the ellipse. 

The general solution of the elliptic two-body problem 
is given in terms of the time t and six integration con- 
stants called the orbital elements. Three of the latter 
are associated with inherent properties of the ellipse: the 
semimajor axis length o, the eccentricity e, and a phase 
variable Iq called the mean anomaly at the epoch. The 
other three describe the position of the ellipse in the ref- 
erence coordinate system: the inclination i, the longitude 
of the ascending node il, and the angle uj defined above. 
The orbital frequency is denoted n, and it is related as 
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FIG. 6: Definitions for the elliptical orbit. 



usual to the period Pi, = 27r/n. Standard expressions 
relate these various quantities. For example, the eccen- 
tric anomaly E is related to the mean anomaly lo at the 
epoch t = to hy 

E -esinE = lo + n{t-to). (161) 

The reader is referred to Refs. [75] for further details. 

2. Secular changes 

We begin the analysis with the derivation of the sec- 
ular changes in the orbital elements of a binary-pulsar 
system that arise from Lorentz violation. Adopting the 
point-mass limit as described above, Eq. (52) can be used 
to determine the relative coordinate acceleration of the 
pulsar and companion, which are taken to be located at 
fi and r2, respectively. The result can be written in the 
form of Eq. (99), with the quadrupole moment set to zero 
and the effects from external bodies omitted. 

It is convenient for the analysis to choose a post- 
newtonian frame with origin at the center of mass of 
the binary-pulsar system. The newtonian definition of 
the center of mass suffices here because the analysis is 
perturbative. The frame is taken to be asymptotically 
inertial as usual, so that the coefficients for Lorentz vi- 
olation obey Eqs. (11). The coordinates in this frame 
are denoted = {t,x^), and the corresponding basis 
is denoted e^. The issue of matching the results to the 
Sun-centered frame is addressed in Sec. VE5 below. 

In this post-newtonian binary-pulsar frame, the rel- 
ative acceleration of the pulsar and companion can be 
written 



2G6m 



where the definitions of r, r, il, M and 5m parallel those 
adopted in Table 1 of Sec. VB. The ellipses here rep- 
resent corrections from multipole moments, from general 
relativistic effects at 0(4), and from Lorentz- violating ef- 
fects at 0(4). Note that enhancements of secular effects 
stemming from general relativistic corrections at 0(4) 
are known to occur under some circumstances [74], but 
the issue of whether 0(4) enhancements can occur in 
the SME context is an open question lying beyond our 
present scope. 

To establish the observational effects of Lorentz viola- 
tion that are implied by Eq. (162), we adopt the method 
of osculating elements [75]. This method assumes the 
relative motion of the binary pulsar and companion is 
instantaneously part of an ellipse. The subsequent mo- 
tion is obtained by letting the orbital elements vary with 
time while keeping the equations of the ellipse fixed. In 
the present instance, the unperturbed ellipse solves the 
equation 



GM 



(1 + |s°°)r^- 



(163) 



Note that this implies the frequency and the semimajor 
axis of the orbit are related by 



GM(1 
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(164) 



The remaining terms are then viewed as a perturbing 
acceleration a'-' , given in this case by 
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s'^^v^r"]. (165) 



+ — 3— [s°'=t;'=r^' - f^v'^r''] +...., 



(162) 



The idea is to extract from this equation the variations 
in the six orbital elements a, e, Iq, i, fl, lo describing the 
instantaneous ellipse for the relative motion. 

The elliptical orbits of the pulsar and companion, re- 
ferred to the chosen post-newtonian frame, are illustrated 
in Fig. 7. The relative motion can also be represented by 
an ellipse in the same plane. To specify the orientation 
of the orbit, three linearly independent unit vectors are 
needed. A natural set includes one vector fc normal to 
the orbital plane, a second vector P pointing from the 
focus to the periastron, and a third vector Q = k x P 
perpendicular to the first two. This triad of vectors can 
be expressed in terms of the spatial basis for the post- 
newtonian frame. Explicitly, these vectors are given by 

P = (cosficosw — coszsinr2sinti;)ei 

-|- (sin n cos u) + cos i cos sin w) 62 

-l-sinisinwea, 
Q = — (cosOsino; -|- cosisinf2cosa;)ei 

-|- (cos i cos ft cos w — sin f2 sin w) 62 

-|-sinicoswe3, 
k = sin z sin Qei — sin i cos ^62 

+ cosies. (166) 
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FIG. 7: Elliptical orbits in the post-newtonian frame of the 
binary pulsar. 



In terms of these vectors, the unperturbed eUiptical orbit 
can be expressed as 



r = 



a(l 



1 + e cos / 



(Pcos/ + gsin/). 



(167) 



Following standard methods, we insert the solution r 
for the elliptic orbit given by Eq. (167) into the expres- 
sion (165) for the perturbative acceleration, project the 
result as appropriate to extract the orbital elements, and 
average over one orbit. The averaging is performed via 
integration over the true anomaly /, which serves as the 
phase variable. Details of this general procedure are dis- 
cussed in Refs. [75]. 

After some calculation, we find the secular evolutions 
of the semimajor axis and the eccentricity are given by 
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where the eccentricity function e{e) is given by 



s{e) = l-(l-e2)V2. 



(168) 



(169) 



The secular evolutions of the three orbital elements as- 
sociated with the orientation of the ellipse relative to the 



post-newtonian frame are 
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(170) 



Finally, the secular evolution of the mean anomaly at the 

epoch is 
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(171) 



In all these expressions, the coefficients for Lorentz vio- 
lation with subscripts P, Q, and fc are understood to be 
appropriate projections of 's'^^ along the unit vectors P, 
Q, and ^, respectively. For example, sp = s^^P^, while 



P'^Q^Q''. 

The analysis of the timing data from binary-pulsar sys- 
tems permits the extraction of information about the sec- 
ular changes in the orbital elements. In the present con- 
text, performing this analysis for a given binary pulsar 
would permit measurements of the projected combina- 
tions of coefficients for Lorentz violation appearing in the 
equations above. This is discussed further in Sec. VE4. 
Matching the projections of the coefficients to standard 
components in the Sun-centered frame can be performed 
with a Lorentz transformation, described in Sec. VE5. 



3. Timing formula 

Accompanying the above secular changes in the or- 
bital elements are modifications to the standard pul- 
sar timing formula arising from Lorentz violation in the 
curved spacetime surrounding the binary pulsar. The ba- 
sic structure of the pulsar timing formula can be derived 
by considering the trajectory of a photon traveling from 
the emitting pulsar to the Earth. Together with the rela- 
tion between proper time at emission and arrival time on 
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Earth, this trajectory is used to determine the number 
of pulses as a function of arrival time. For simplicity, the 
derivation that follows neglects any effects on the photon 
trajectory from gravitational influences of bodies in the 
solar system and effects from the motion of the detector 
during measurement. 

The photon trajectory is determined by the null condi- 
tion ds^ = 0. The relevant metric is the post-newtonian 
result of Eqs. (35), (36), (37) taken at 0(2) and incor- 
porating effects of the perturbing companion of mass m2 
at r2- In practice, the time-delay effects from the pulsar 
itself produce only a constant contribution to the result 
and can be disregarded. The null condition can then be 
solved for dt and integrated from the emission time tem 
at the pulsar located at fi to the arrival time tarr at the 
Earth at fs- Assuming \rE\ ^ |ri| and \ fE\ ^ \f2\, we 
find the time delay farr — iem for the photon is given by 
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(172) 



This result is a simplification of the general time-delay 
formula, which is derived in Sec. VF2 below. In Eq. 
(172), te = I^eI is the coordinate distance to the Earth, 
r = ri — r2 is the pulsar-companion separation as before, 
and n is a unit vector along fs — ri . The quantity A-'^ 
contains combinations of coefficients for Lorcntz violation 
and n, given explicitly in Eq. (198) of Sec. VP 2. 

We seek an expression relating the photon arrival time 
to the proper emission time. The differential proper time 
dr as measured by an ideal clock on the surface of the 
pulsar is given in terms of the differential coordinate time 
dt by the expression 



dr = dt[l -{1 + 
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The potentials appearing in this expression include those 
from both the pulsar and its companion, and they are 
to be evaluated along the worldline of the ideal clock. 
It suffices here to display the explicit potentials for the 
companion, so that 
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where Ui and are the contributions from the pulsar. 
In Eq. (173), v is the velocity of the clock, which can 
be identified with the velocity vi of the pulsar. Using 
the time derivative of Eq. (167) to obtain vi and explic- 
itly evaluating Eq. (173) along the pulsar trajectory, we 
obtain 
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All effects from the pulsar itself contribute to constant 
unobservable shifts and have been neglected in this result. 
Note that the factors of (1 -|- |s°°) appearing in Eq. (175) 
merely act to scale G and hence are unobservable in the 
present context. For simplicity in what follows, we absorb 
them into the definition of G. 

To integrate Eq. (175), it is convenient to change vari- 
ables from the coordinate time to the eccentric anomaly 
E. This makes it possible to insert the expression 



r = a(cos E — e)P + a(l 
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into Eq. (175) and perform the integration over E. Up 
to constants, we obtain 
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The functions J^i and J^2 appearing in the last two terms 
are given by 
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Note that in the circular limit e — > the result (177) 

remains finite, despite the seeming appearance of poles 
in e. The constant in the definition (161) is fixed by 
requiring 



E — esinE = nitun — te)- 
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The integration assumes that the positions of the pulsar 
and its companion are constant over the photon travel 
time near the binary-pulsar system. 

The desired expression for the proper time t is ob- 
tained by combining the result (172) with Eq. (177) to 
yield an expression as a function of the arrival time. It 
turns out that the Lorentz-violating corrections to the 
Shapiro time delay, which arise from the last term of Eq. 
(172), are negligibly small for the cases of interest and so 
can be disregarded. We find 

T = te^rr - A{cosE - e) - {B + C)smE 
27T (AsinE- BcosE) 
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In this expression, ai = am^/M and the quantities A, B, 
and C are given by 



A = 
B = 

C = 



ai sm « smo;, 
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The reader is reminded that all these equations are to 
be used in conjunction with the method of osculating 
elements, for which the orbital elements vary in time. 

In practice, the application of the above results in- 
volves inserting the result (180) into an appropriate 
model for the number N{t) of emitted pulses as a func- 
tion of pulsar proper time r. The model for N(t) is 
usually of the form 



N(t) = NQ^-yT + + . . 



(182) 



where v is the frequency of the emitted pulses. The re- 
sult of the substitution is an expression representing the 
desired timing formula for the number of pulses N{ta,ir) 
as a function of arrival time. 



4- Experiment 

We are now in a position to consider some experimental 
implications of our results. Consider first the timing for- 
mula (182) with the substitution of the expression (180) 
for the proper time. Modifications to this formula aris- 
ing from Lorentz violation can be traced to two sources. 
First, the terms in Eqs. (180) and (181) proportional to 
spp, sqq, and spq represent direct corrections to the 
proper time arising from Lorentz violation. Second, the 
structural orbital elements entering the expression for the 
proper time r acquire secular Lorentz-violating correc- 
tions, given in Eqs. (168) and (170). Both these sources 
of Lorentz violation affect the predictions of the timing 
formula. 

In applications, the timing formula N can be regarded 
as a function of a set of parameters 



N = N{Pb,Pb,e,e,w,(j,...). 



(183) 



One approach to measuring the coefficients for Lorentz 

violation would be to fit the observed number of pulses 
A'^ as a function of arrival time iarr using a least-squares 
method. In principle, an estimate of the sensitivity of 
this approach could be performed along the lines of the 
analysis for the conventional timing formula in Ref. [76] . 
However, we adopt here a simpler method that nonethe- 
less provides estimates of the sensitivities that might be 
attainable in a detailed fitting procedure. 

To obtain these simple estimates, wc limit attention to 
the secular changes in the orbital elements e and u> given 
in Eqs. (168) and (170). Values for e and di obtained 
from the timing analysis can be used to probe combina- 
tions of coefficients in s'^'^ . Assuming only measurements 



from the variations of the orbital elements are considered, 
the number of independent coefficients for Lorentz vio- 
lation to which binary-pulsar systems are sensitive can 
be quantified. Through the measurement of w and e, a 
given binary-pulsar system is sensitive to at least 2 com- 
binations of coeflacients. These are 
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In principle, sensitivities to the secular changes in i and 
Q, might increase the number of independent coefficients 
that could be measured in a given system. 

Since there are 9 independent coefficients in s'^" , it 
is worth studying more than one binary-pulsar system. 
Typically, distinct binary-pulsar systems are oriented dif- 
ferently in space. This means the basis vectors P, Q, 
k differ for each system, and hence the dependence on 
the coefficients for Lorentz violation in the Sun-centered 
frame also differs. Examination of Eq. (184) reveals that 
the combination spp -\- sqq is inaccessible to any orien- 
tation and that 5 independently oriented binary-pulsar 
systems can in principle measure 8 independent coeffi- 
cients in s^^ . If information from i and 51 can also be 
extracted, then fewer than 5 binary-pulsar systems might 
suffice. 

The error bars from existing timing analyses offer 
a rough guide to the sensitivities of binary pulsars to 
Lorentz violation. As an explicit example, consider the 
binary pulsar PSR 1913 -|- 16. The errors on e and 6j are 
roughly 10^-'^^ and 10~^, respectively [77]. These imply 
the following sensitivities to Lorentz violation: Sg < 10~^ 
and Suj < 10"-^^. Interesting sensitivities could also be 
achieved by analysis of data from other pulsar systems 
[78, 79]. 



5. Transformation to Sun-centered frame 

In the preceding discussions of tests with binary pul- 
sars, all the coefficients for Lorentz violation enter as pro- 
jections along the triad of vectors P, Q, k defining the 
elliptical orbit. These projected coefficients involve the 
post-newtonian frame of the binary-pulsar system. The 
explicit construction of the transformation that relates 
these projections to the standard coefficients in the 
Sun-centered frame is given next. 

The first step is to express the various projections of 
the coefficients directly in terms of the orbital angles i, 
and oj and the coordinate system x'^ of the binary pulsar. 
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This involves substitution for the projection vectors using 
Eqs. (166). For example, for the projection spp we have 

Spp = P^P^s^^ 

= s^^(cosOcosa; — cosisinnsino;)^ 
+s^^(sinncosa; + cosicosf2sina;)^ 

+s^^ sin^ i sin^ w 

+2s^'^ (cos 51 cos ui — cos i sin ft sin iv) 

X (sin Q cos ui + cos i cos sin w) 
+2s^^ sin i sin u> (cos Q cos w — cos z sin O sin w) 
+2s^^ sin i sin w (sin O cos u + cos i cos f2 sin a; ) . 

(185) 

The resulting expressions give the projected coefScients 
in terms of the coefficients s^'^, with components in the 
basis of the binary pulsar. 

The conversion to coefficients in the Sun-centered 
frame can be accomplished with a Lorentz transforma- 
tion. We have 

s"^ = A^^A V^n. (186) 

The Lorentz transformation consists of a rotation to align 

the axes of the binary-pulsar and Sun-centered frames, 
along with a boost to correct for the relative velocity. The 
latter is taken as approximately constant over the time 
scales of relevance for the experimental observations. 

The components of the Lorentz transformation are 
given by 

A°T = 1. 

hi J = RjJ. (187) 

In this expression, 0'^ is the boost of the binary-pulsar 
system as measured in the Sun-centered frame. The ma- 
trix with components i?^"^ is the inverse of the rotation 
taking the basis for the Sun-centered frame to the 
basis e^^ for the binary-pulsar frame. Explicitly, these 
components are 
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In these expressions, the angles a, (3, and 7 are Euler 
angles defined for the rotation of ej to as follows: 

first, rotate about ez by 7; next, rotate about the new 
Y axis by a; finally, rotate about the new Z axis by 7. 



F. Classic Tests 

Among the classic tests of general relativity are the 
precession of the perihelia of the inner planets, the time- 
delay effect, and the bending of light around the Sun. 
The sensitivities to Lorentz violation attainable in ex- 
isting versions of these tests are typically somewhat less 
than those discussed in the previous sections. Nonethe- 
less, an analysis of the possibilities reveals certain sen- 
sitivities attaining interesting levels, as described in this 
subsection. 



1. Perihelion shift 

To obtain an estimate of the sensitivity to Lorentz vi- 
olation offered by studies of perihelion shifts, the results 
from Sec. VE can be adapted directly. The reference 
plane can be chosen as the ecliptic, so the inclination i is 
small for the inner planets. 

We seek the perihelion shift with respect to the 
equinox. The angular location uj of the perihelion with 
respect to the equinox is given hy uj = uj + cos i. The 
change in lu per orbit can be obtained by appropriate use 
of Eqs. (170) with dm/M w 1. Neglecting terms propor- 
tional to sini, we find that the shift in perihelion angle 
per orbit is 

„ r(e^-2£),_ _ , 2na(e2-e)_i 

(189) 

Substituting the values of e and na for Mercury and 
the Earth in turn, we find the corresponding perihelion 
shifts in arcseconds per century (C) are given by 

ci)^ « 7 X W'^"s^C-\ 

(i)® » 2 X 10^"s®C-^ (190) 

In these expressions, the perihelion combinations of SME 
coefficients for Mercury and the Earth are defined by 

« {spp - Sqq) - 6 X 10"^SQ, 

s© w {spp-SQo) - 5 X IQ-^SQ. (191) 

The reader is cautioned that the vector projections of 
along P and Q differ for Mercury and the Earth. 
Indeed, Eqs. (166) reveal that they differ by the values of 
i, LO and f2, once again indicating the value of considering 
more than one system for measurements separating the 
independent values of s''''. 

The error bars in experimental data for the perihelion 
shifts of Mercury and the Earth provide an estimated 
attainable sensitivity to Lorentz violation. For example, 
adopting the values from Refs. [30, 72] , the error bars on 
the perihelion shift of Mercury and the Earth are roughly 
given by 0.043" and 0.4" C"\ respectively. Taking 
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these values as an upper bound on the size of the observ- 
able shifts in Eq. (190), we obtain attainable sensitivities 
of < 10-9 and s© < 10"^ 

For simplicity in the above, attention has been focused 
on Mercury and the Earth. However, the varying ori- 
entations i, u, O of other bodies in the solar system 
suggests examining the corresponding perihelion shifts 
might yield interesting sensitivities on independent coef- 
ficients for Lorentz violation. For example, observations 
of the asteroid Icarus have been used to verify to within 
20 percent the general-relativistic prediction of 10" C"^ 
for the perihelion shift [80]. Translated into measure- 
ments of Lorentz- violating effects, this would correspond 
to a sensitivity at the 10~^ level for the combination 

Sic = (SPP - Sqq) - 7 X 10~*SQ. 



2. Time-delay effect 

For the binary pulsar, the Lorentz- violating corrections 
to the Shapiro time-delay effect are negligible compared 
to other effects, as described in Sec. VE. However, solar- 
system tests involving the time-delay effect have the po- 
tential to yield interesting sensitivities. To demonstrate 
this, we derive next the generalization of the usual time- 
delay result in general relativity to allow for Lorentz vi- 
olation. 

We seek an expression for the time delay of a light 
signal as it passes from an emitter to a receiver in the 
presence of a massive body such as the Sun. Assuming 
that light travels along a null geodesic in spacetime, the 
photon trajectory satisfies the condition ds^ = 0. In 
terms of the underlying post-newtonian coordinates, this 
can be written in the form 



~dt dTj 



1/2 



= 1 



00 



dx^ dx^ 
~dt~dr' 

(192) 



The metric components are taken to be the 0{2) forms 
of Eqs. (35), (36), and (37) in the presence of Lorentz 
violation and with a perturbing point mass m2 located 
at r2. 

Solving Eq. (192) for dt, we integrate from an initial 
emission spacetime point (tcm, n.) to a final arrival space- 
time point (tarrjre)- At order 0(2), the zeroth-order 
straight-line solution 



r{t) = n{t - Um) + ri(iem) 



(193) 



c;an be inserted into the metric components in Eq. (192). 
Here, n is a unit vector in the direction te — ri. In 
performing the integration, the positions of the emitter 
at r*!, receiver at rg;, and perturbing body at ^2 a-r^^ 
taken to be constant during the travel time of the light 
pulse. 

The resulting expression for the time difference farr — 
iem is found to be 

tarr - ^em = |rE - fi | Ai A2(iarr) " A2(iem)- (194) 



The terms Ai and A2 on the right-hand side are given 

by 



Ai = 



l^^arr ~t- C(^arr)| ^ ^arr ^ ^ * C(^arr) 



X In 



I'^^em + C(iem)| + t 

em + n ■ (item) 



(195) 



and 

A2(t) 



Gm2A^'= 



h^n^iiin-CY -d')t + n-a'] 



d2|nt + C| 

- 2We{^ + h . Ct) + CCHt + n-0], 

(196) 

where ( and d are functions of a time argument and are 
defined as 

C(*) = n(iem) - r2{t) - ntem, 

d{t) = n X (C X n). (197) 

The combination of coefficients A-''^, which includes the 
effects due to Lorentz violation, is given by 

^jfe = 2d^'' + Ai'', (198) 



where 



if = (sOO-sO'^^i)5ife 



(199) 



Note that the trace part of A-''^ cancels in Eq. (196). 

The result (194) is the generalization to the pure- 
gravity minimal SME of the standard time-delay formula 
of general relativity, which can be recovered as the limit 
-gtj.v _ Q principle, Eq. (194) could be applied to solar- 
system experiments including, for example, reflection of 
a signal from a planet or spacecraft. A detailed analysis 
along these lines would be of definite interest but lies be- 
yond our present scope. However, the form of Eq. (194) 
indicates that differently oriented experiments would be 
sensitive to different coefficients for Lorentz violation. 

At present, the most sensitive existing data on the 
time-delay effect within the solar system are obtained 
from tracking the Cassini spacecraft [81]. Assuming a 
detailed analysis of these data, we estimate that sensi- 
tivities to various combinations of the coefficients s^^ 
could attain at least the 10~^ level. This is compara- 
ble to the sensitivities of GPB discussed in Sec. VD3, 
but the differences between the two spacecraft trajecto- 
ries makes it likely that the attainable sensitivities would 
involve different combinations of the coefficients s'^^. 

We remark in passing that a generalization of the stan- 
dard formula for the bending of light should also exist, 
following a derivation along lines similar to those lead- 
ing to Eq. (194). Current measurements of the deflection 



39 



of light by the Sun have attained the 10~^ level [72], 
which suggests sensitivities to Lorcntz violation at this 
level might already be attainable with a complete data 
analysis. Future missions such as GAIA [82] and LATOR 
[83] may achieve sensitivities to coefficients for Lorentz 
violation at the level of 10~^ and 10~*, respectively. 

VI. SUMMARY 

In this work, we investigated the gravitational sector of 
the minimal SME in the Ricmann-spacctime limit, with 
specific focus on the associated post-newtonian physics. 
The relevant basics for this sector of the SME are re- 
viewed in Sec. II A. There are 20 independent back- 
ground coefiicient fields u, s^" , and t'^^t^'^ ^ with corre- 
sponding vacuum values u, s^", and t'^^^" _ 

The primary theoretical challenge to the study of post- 
newtonian physics is the extraction of the linearized field 
equations for the metric while accounting for the dynam- 
ics of the coefficients for Lorentz violation. At leading 
order, this challenge is met in Sec. II B. Under a few mild 
conditions, the linearized effective field equations appli- 
cable for any theory with leading-order Lorentz violation 
are obtained in the practical form (23) . 

Section III is devoted to the post-newtonian metric. 
Its derivation from the effective linearized field equations 
is described in Sec. Ill A. Under the simplifying assump- 
tion u = 0, the result is given in Eqs. (35)- (37). Section 
IIIB obtains the equations of motion (42) and (43) de- 
scribing the corresponding post-newtonian behavior of a 
perfect fluid. The results are applied to the gravitat- 
ing many-body system, with the point-mass equations of 
motion and lagrangian given by Eqs. (52) and (54). The 
post-newtonian metric obtained here is compared and 
contrasted with other existing post-newtonian metrics in 
Sec. Ill C. Finally, to complete the theoretical discussion, 
illustrative examples of the practical application of our 
methodologies and results to specific theories are given 
in Sec. IV in the context of various bumblebee models. 

The largest part of the paper, Sec. V, is devoted to 
exploring the implications of these theoretical results for 
a variety of existing and proposed gravitational exper- 
iments. It begins with a generic discussion of coordi- 
nate frames and existing bounds in Sec. VA. Section 
VB considers laser ranging both to the Moon and to 
artificial satellites, which offers promising sensitivity to 
certain types of Lorentz violation controlled by the co- 
efficients s"^ in the Sun-centered frame. The key result 
for the relative acceleration of the Earth-satellite system 
is given in Eq. (99) and the associated definitions. The 
technical details of some of the derivations in this section 
are relegated to Appendix A. 

Section VC considers Earth-based laboratory exper- 
iments. The expression for the local terrestrial accel- 
eration as modified by Lorentz violation is obtained as 
Eq. (126) in Sec. VCl. Applications of this result to 
gravimeter measurements are treated in Sec. VC 2, while 



implications for experiments with torsion pendula arc de- 
scribed in Sec. VC3. The latter include some proposed 
high-sensitivity tests to the coefficients based on the 
local anisotropy of gravity. 

Possible sensitivities to Lorentz violation in measure- 
ments of the spin precession of orbiting gyroscopes are 
considered in Sec. VD. The generalized formula describ- 
ing the spin precession is found in Eq. (143), and the 
secular effects are isolated in Eq. (153). The results are 
applied to the Gravity Probe B experiment [70]. This ex- 
periment is shown to have potential sensitivity to effects 
involving spin precession in a direction orthogonal to the 
usual geodetic and frame-dragging precessions predicted 
in general relativity. 

Timing studies of signals from binary pulsars offer fur- 
ther opportunities to measure certain combinations of 
coefficients for Lorentz violation. Section V E derives the 
secular evolution of the six orbital elements of the binary- 
pulsar motion and calculates the modifications to the pul- 
sar timing formula arising from Lorentz violation. Fits to 
observational data for several binary-pulsar systems have 
the potential to achieve high sensitivity to coefficients of 
the type s~^. 

Section VF discusses the attainable sensitivities in 
some of the classic tests of general relativity, involving the 
perihelion precession, time-delay, and bending of light. 
Perihelion shifts for the Earth, Mercury, and other bodies 
are considered in Sec. VF 1, while versions of the time- 
delay and light-bending tests are studied in Sec. VF2. 

Table 6 collects the estimated attainable sensitivities 
for all the experiments considered in Sec. V. Each hori- 
zontal line contains estimated sensitivities to a particular 
combination of coefficients for Lorentz violation relevant 
for a particular class of experiment. Most experiments 
are likely to have some sensitivity to Lorentz violation, 
but the table includes only the better estimated sensitiv- 
ities for each type of experiment. 

Relatively few sensitivities at these levels have been 
achieved to date. In Table 6, estimated sensitivities that 
might be attainable in principle but that as yet remain 
to be established are shown in brackets. Values without 
brackets represent estimates of bounds that the theoreti- 
cal treatment in the text suggests are already implied by 
existing data. All estimates are comparatively crude and 
should probably be taken to be valid only within about 
an order of magnitude. They are based on existing tech- 
niques or experiments, and some planned or near-future 
experiments can be expected to supersede them. 

In conclusion, both astrophysical observations and ex- 
periments on the Earth and in space have the potential 
to offer high sensitivity to signals for Lorcntz violation of 
purely gravitational origin. A variety of tests is required 
to span the predicted coefficient space for the dominant 
Lorentz-violating effects. These tests are ciurently fea- 
sible, and they represent a promising direction to follow 
in the ongoing search for experimental signals of Lorentz 
violation from the Planck scale. 
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in the text. The second and third parts are devoted to 
forced and unforced oscillations, respectively. The final 
part extracts the radial oscillations needed for the anal- 
ysis in Sec. VB. 



APPENDIX A: SATELLITE OSCILLATIONS 

This appendix is devoted to calculating the Lorentz- 
violating corrections to the observable relative Earth- 
satellite distance r = [ri — in the Sim-centered frame. 
These corrections arc needed for the analysis of lunar 
and satellite laser ranging in Sec. VB. For simplicity, 
the only perturbing body is taken to be the Sun. Also, 
we limit attention to the dominant frequencies of direct 
relevance for the text. A more detailed analysis could 
reveal other frequencies that would permit sensitivity to 
different combinations of coefficients for Lorentz viola- 
tion. 

The derivation is separated into four parts. The first 
discusses generalities of the perturbative expansion and 

establishes notational conventions. Some of the results 
in this part are also applied in other contexts elsewhere 



1. Generalities 

The procedure adopted here involves perturbation 
techniques similar to those presented in Refs. [52, 84]. 
The equation (99) for the coordinate acceleration of the 
Earth-satellite distance is expanded in a vector Taylor 
series about an unperturbed circular-orbit solution Tq . 
This unperturbed circular orbit satisfies 



^0 =aoK), 



(Al) 



where 



with 



, GMr-^ nlr-^ ZQr-' 



GMq 



(A3) 
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The frequency of the circular orbit is defined as 



where the oscillating basis vectors Pv, f,^, are given by 



GM 



2"e 



3Q 
2r5- 



(A4) 



Contributions from perturbing bodies and the Earth's 
quadrupolc potential are included to mute resonant ef- 
fects [48, 51, 52], arising in this instance from Lorentz- 
violating accelerations. 

Denoting the perturbation by e'', we find 

-{e''dK)a'o{ro) - ^{e^e'^dKdL)a-^{ro) + . . . 

+ai,{ro) + (e-^aif)a^,(ro) + . . . 
+a^,(ro) + (e^5K)a^.(ro) + ... 

+aiv{ro) + {e^dK)aiy{ro) + .... 

(A5) 

The ellipses hero represent higher-order terms in the vec- 
tor Taylor series. The quantities a'j., and a'q, are given 
by modifications of Eqs. (101) and (102), respectively, in 
which the tidal and quadrupole terms for Eq. (A2) are 
removed. The term aly(r) is defined in Eq. (103). 

It is convenient to decompose Eq. (A5) using a basis of 
3 ortlionormal unit vectors p, f , &. These vectors arc, re- 
spectively, normal, tangential, and perpendicular to the 
unperturbed orbit. In the basis for the Sun-centered 
frame, the explicit form for p is 

p = [ cos a cos{coT + 6) — sin a cos /3 sin(a;T + 9) , 
sin a cos{ijT + 9) + cos a cos /3 sin(a;T + 9) , 
sin/Jsin(wT-|-6i)), (A6) 

where 9 is the phase of the orbit with respect to the line 
of ascending nodes. Explicit forms for the other two basis 
vectors can be obtained from the equations 



and 



LOT — dp/dT 



a = p X T. 



(A7) 



(A8) 



Note also that rg = rop^ , where vq is the radius of the 
circular orbit. 

For our purposes, it suffices to assume that m2 » mi 
and to take R'^ as the vector pointing to the Earth's 
location. The Earth's motion can also be approximated 
as a circular orbit around the Sun. This means that 

R = ( — cosJIqT, — cos77sinO0T, — sinr/sinfi^T) 

(A9) 

in the Sun-centered frame, so we can write R^ = RR^ 
and = dR-'/dT. 

We write the perturbative oscillations in the general 
form 



Pi, = pcos{vT + 9v), 
TV = Tsm{i>T + 9^), 
Gv = acos{uT + 9i,). 



(All) 




With respect to this oscillating basis, Eq. (A5) can be 
expressed as a matrix equation. For example, takes 
the column-vector form: 



(A12) 



The types of oscillations contained in Eq. (A5) can be 
split into two categories: forced or inhomogcneous oscilla- 
tions, and unforced or homogeneous oscillations [52, 84]. 
The forced oscillations arise primarily from the terms 
ax', dq', cthv- The unforced oscillations are the nat- 
\iral oscillations of the Earth-satellite system in the pres- 
ence of the perturbative forces. They are controlled by 
the matrix structure of the derivative terms in Eq. (A5), 
such as e^e^^K^iao (?'o)/2. Next, we study each of these 
categories of oscillations in turn. 



Forced oscillations 



The forced oscillations due to Lorentz violation em- 
anate primarily from the terms in cSlv. For oscillations 
near the resonant frequency u), however, the tidal and 
quadrupole derivative terms in Eq. (A5) can play a sig- 
nificant role. Keeping all the relevant terms, the lowest- 
order contributions arise from the equation 



(e^5K)ai((ro) = ai^v(ro) + (e^9K)a^-(ro) 

+ {e''dK)ai^,{ro). (A13) 



(AlO) 



In the oscillating basis of Eq. (A12), the relevant part of 
this equation takes the matrix form 

K (y) ■ e^fi^ = (ax + "q) • el-.e, + diM{v, 9^) 

(AM) 

for each frequency v. Here, we define K (i^) to be the ma- 
trix operator with components 5 {d^ / dT"^) — i9ifao (ro) 
in the Sun-centered frame, while ax and uq denote the 
diagonal parts of Okcxt' ^^'^ ^kchq,, respectively. 

By evaluating all terms in Eq. (A13) in the oscillat- 
ing basis, we find that the components of these matrix 
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operators for a given frequency v are given by 



K{v) = - 



( + 3a;2 + 3^2^ 2ujv 

,2 



V 



2ujv 




V 

_ ^2 



30a 




/ ai - i ^ 

ai - i 

V a2-i/ 

21sin2/? 










3 c-in-^ , 



-3-|sin^/3 y 



(A15) 



Here, the quantity O^ is defined by 



OS 



O^ +^ 



(A16) 



by the expressions 

aLv(2w,26'- 0-20,) = 

aLv(w,6' + cra,,i) = 

aLv(w,^ - C7(^,2) = 



aLv(0©,crnffi,2) = 




-GMV^ I I . (A19) 



2r2 



while ai and a2 are given by 

ai = \ [cos^ Q!(l + cos^ 13 cos^ ??) + sin^ j3 siv? rj 

+ sin^ Q!(cos^ (3 + cos^ 77) 

+2 cos a cos /3 sin (3 cos 77 sin 77] , 
0-2 = \ [sir? a s\v? [3 + cos^ a sin^ (3 cos^ 77 

+ cos^ (3 sin^ 77) — cos a cos /3 sin (3 cos 77 sin 77. 

(A17) 



Inspection of the form of q:lv given in Eq. (103) re- 
veals that the dominant oscillations forced by the Lorentz 
violation occur at the frequencies 2w, w, and O©. The 
amplitudes at these frequencies are determined by the 
following matrix equations: 



K (2w) 

K (Oe) • ehgj,<^„ 



(ax + aq) 
+aLv(2w,26'-(T2,.), 

(ax + aq) ■ eL,e+<7„,i 
+aLv(w,^-|-crt^,i), 

(ax + aq) ■ eL,e-<T„,2 
+aLv(i^- - 0-^^,2), 

(ax + aq) • ebg,,^^^,! 
+aLv(0©, aO(B,i), 

(ax + aq) • ene.^siffi,2 
+QLv(0®,ane,2). (A18) 



The various combinations of the coefficients for Lorentz 
violation have magnitudes and phases defined by 



S2uj COS a2ui 

S20J sin a2ui 

Su,l C0SC7a,,l 

Su,i sincTi^^i 

Sw,2 C0S(Ta,,2 

Sa,,2sinc7a,,2 
sne,i cos crn^,! 
sn®,! sin an^^i 
sn®,2 cosc7ne,2 



i(^"-s-), 



^02 



= s 



23 



'T'V . 

COS jjs + sm rjs , 

sin a cos 771'^^ 



(sin /3 sin 7; + cos Of cos /? cos a) s , 

(A20) 



^02 



su^,2sma^^^2 = cos as -sin a cos /3s 



In these expressions, the quantities s^'^ are combinations 
of the coefficients for Lorentz violation in the Sun- 
centered frame. For s^^ — s^^, s^^ , s"^, and s*^^, these 
combinations can be found in Eqs. (107) of Sec. VB, 
while for s^^ and s^^ they are given by 



isin/3sin2a(s^-^ 
+ cos /3(cos as"^^ -f 



s ) — sin /3 cos 2as 



XY 



sm as 



YZ 



s = 



-isin2/3[i(s^^ + s^^)-s^^] 
-1-5 sin 2/3 cos 2a(s^^ -s 



+i sin 2/3 sin 2as 



— sin a cos 2/3s -|- cos a cos 2/3s 



(A21) 



In these equations, the relevant terms in iSlv are given 



Using the appropriate inverse matrix, the solutions for 
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2a; and can be written as 



Unforced oscillations 



e2w,6l-(72u< 



S2Loro 



en 




(A22) 



where 



bi = i - 3ai + 



21Qsin2/3 



+9ai 



630 sin^ p 



(A23) 



with fli defined in Eq. (A17). 

Special consideration is required for the near-resonant 
frequency u). The matrix to be inverted for this case is 
K'{lu) = K{lo) — q:t — q:q- To sufficient approximation, 
the determinant of this matrix can be expanded in terms 
of the anomalistic frequencies ojq and Wq. These are as- 
sociated with eccentric motions in the plane of the orbit 
and perpendicular to the orbit, respectively. The motion 
involving the phase is confined to the plane of the 

orbit and hence is at frequency wq, while that involving 
the phase 6 — (7^1,2 is perpendicular to the orbit and is at 
frequency lOq. We make use of the approximation 



detK'{uj) = detK'{LJo+uj-iOo) 
w det K'{u)o) + {w — u)o) 



ddctK'iuio) 



(A24) 

along with the similar relation for the frequency w'q. This 
permits the inverse of K'{w) to be determined as 



[K'{uj)\- 



2ijJq{ijJ - Wo) 



/-I 

2 

V 



2 

-4 




where 



7 



uJo{ijj - Wo) 



In consequence, the w amplitudes now read 



(A25) 



(A26) 



uj 5mSu,iVoro 
{(J - Wo) M 



(w - uj'o) 2 



( 



V-1 



(A27) 



Among the unforced oscillations satisfying Eq. (A5) are 
the natural eccentric oscillations of the systc;ni occurring 
in the absence of Lorentz violation. These satisfy a set of 
coupled matrix equations that can be derived from Eq. 
(A5). Here, we are interested in unforced oscillations at 
the anomalistic frequency luq- 

In a simplified scenario, where contributions from the 
right-hand side of Eq. (A5) can be ignored, the anomal- 
istic oscillation satisfies 



K (wn) 



j) 



(A28) 



Requiring the vanishing of the determinant of the matrix 
K' = K — ar — aQ allows the extraction of approximate 
solutions for wq. This gives three solutions, two of which 
are physical. One is denoted coq and corresponds to ec- 
centric motion in the plane of the orbit. The other is 
denoted lu'q and corresponds to eccentric motion perpen- 
dicular to the orbit. Explicitly, we find the approximate 
solutions 



W — Wo 



W — Wn 



2w' 
2w2 



(A29) 



where 



nl = 0|[15ai -6] -I- 



3Q 



0|[3a2 



3Q 



(l + isin2/3). (A30) 



More accurate solutions require a detailed study of a 
system of coupled equations involving other frequencies. 
For example, the eccentric oscillations in the lunar case 
are at frequencies wo and 2(w — f2) — wo [52]. Also, in the 
case of a generic satellite orbit, the Earth's quadrupole 
terms can be expected to play a significant role. 

For our purposes, the primary interest is in the eigen- 
vector for Wo, which is given by 



ero 




(A31) 



The overall normalization of this eigenvector is defined 
in terms of the eccentricity of the orbital motion. 

Using Eq. (A31), the leading unforced oscillation due 
to Lorentz violation can be established. Typically, the 
dominant terms occur at frequencies near the resonant 
frequency, which here is the circular orbit frequency w. 
In this case, we find that the oscillation at frequency 2w — 
Wo is particularly enhanced. In the Sun-centered frame, 
the dominant terms contributing to the amplitude are 
the frequency-mixing terms in e^e^dKdLCtQ{ro)/2 and 
e^dKCtiM that lower the frequency by 2w. The former is 
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obtained by substitution of the forced-oscillation solution 
for ^2uj e-cr^ (^22), while the latter is obtained 

from 



2u\K 
J 



(A32) 



2r5 



The resulting matrix equation takes the form 



K (2a; — wq) • e2u-uo,29-<72u, = 

(A33) 

where the operators on the right-hand side are defined in 
the oscillating basis by 



2 a, 



-2w S2w^ 



-2a! 



\ 



1 3 

/i 1 
1 I 
Vo -I 



(A34) 



(A35) 



Since the frequency 2a; — ujq is near resonant, the pro- 
cedure to find the inverse of K'{2u) — lvq) can be modeled 
on that used for the frequency u in Eqs. (A24)-(A25). 
We find 



Radial oscillations 



At this stage, wc can extract the information appearing 
in Table 2 for radial oscillations at the frequencies lo, fi®, 
2a;, 2u) — u>o. The procedure is to take the p component 
of each column vector ?u.0,,, to multiply by cos{vT + O^,)^ 
and to expand the result in terms of sines and cosines 
using the phase definitions in Eq. (A20). 

After some calculation, we find the following radial os- 
cillations: 



{e2o,)p = -j^(s" -s'')roCOs(2a)T + 20) 
-is^^ro sin(2a;T-^26'), 

/ N W (Sm)vnrn .no . „ „s 

{ejp = , ^''Z ° [s°^cos(a;r + ^) 



(a; - a;o) M 



-s"' sin(a;T-F6i)], 



(en®)p = — VEro[snc cos OeT-FsQ^ sin figjT], 

02 



(UJ - LtJo) 
' 1 -22 



(i(s" - s22) cos[(2a; - uJo)T + 29] 
-hs^^sin[(2a;-wo)r + 26']). (A38) 



[K'{2uj-iOo)r^ 




(A36) 



where 7 is defined in Eq. (A26). It then follows that the 
solution for the oscillations at frequency 2a; — ivq is given 

by 



e2w-wo, 20-0-2^ — 



S2u^ero 
A{w - Wo) 



2 

1 

V 



(A37) 



The new combinations of coefficients for Lorcntz viola- 
tion appearing in these equations are defined by 



-3sn(g,i costTn^,! — ^sn^^2 cosaa^^2 



-SsQ^^isinaa^^i — 550^,2 sin ao^g, 2 



(A39) 
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